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The  Small  Dispersion  Limit    of    the  Korteweg-deVries  Equations 
Peter  D.   Lax  and  David  Levermore 

Abstract 


The  scattering  transform  method  is  used  to  study  the  weak  limit  of 
solutions  to  the  initial  value  problem  for  the  Korteweg-deVries  (KdV) 
equation  as  the  dispersion  tends  to  zero.  In  that  limit  the  associated 
Schrttdinger  operator  becomes  seralclassical ,  so  the  exact  scattering 
data      is     replaced      by      its      corresponding        WKB        expressions.  Only 

nonpositive  initial  data  are  considered;  in  that  case  the  limiting 
reflection  coefficient  vanishes.  The  explicit  solution  of  Kay  and 
Moses  for  the  reflect ionless  inverse  transform  is  then  analyzed,  and 
the  weak  limit,  valid  for  all  time,  is  characterized  by  a  quadratic 
minimum  problem  with  constraints.  With  the  aid  of  function  theoretical 
methods,  solutions  of  this  minimum  problem  are  constructed  in  terms  of 
solutions    to   an   initial  value   problem   for   some   auxiliary  functions. 

The  weak  limit  satisfies  the  KdV  equation  with  the  dispersive  term 
dropped  until  the  finite  time  at  which  its  derivatives  become  infinite. 
Up  to  that  time  the  weak  limit  is  a  strong  L^  limit.  At  later  times 
the  weak  limit  is  locally  described  by  Whitham's  averaged  equations  or, 
more  generally,  by  the  equations  found  by  Flaschka  et  al.  using 
nultiphase  averaging.  For  large  times,  behavior  of  the  weak  limit  is 
studied   directly  from   the    minimum  problem. 
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0.    The   Zero  Dispersion  Limit    for   the  KdV   Equation 

Introduction.      In   this   study  we  analyze   the   behavior  of      solutions 
u(x,t;e)   of 


(1)  u,    -   6uu^  +  e2u^^^=   0 


as  e   ->•   0  while   the   initial  values   are    fixed: 
(2)  u(x,0;e)    =  u(x)    . 

It  is  known  from  computer  studies  that  for  t  greater  than  a  critical 
time,  independent  of  e,  dependent  only  on  the  initial  data,  u(x,t;e) 
becomes  oscillatory  as  e  -»•  0.  The  wave  length  of  these  oscillations  is 
0(e),    their  amplitude    independent   of  e.      This    shows    that 


(3)  lira  u(x, t ;e ) 

e+0 


exists   only    in   the  weak,  sense. 

Our  strategy  for  studying  the  limit  (3)  is  similar  to  one  employed 
by  Hopf  [5]  and  Cole  [1]  to  study  a  similar  but  simpler  singular 
perturbation.   Burgers'   equation: 


(4)  ■  Uj.   +  uu^  =  e   u^^   ,      e    >   0  . 


Hopf    and  Cole   introduced  a  new  variable  V,    defined  by 


(5)  U^  =  u    ,        U    =  -   2e    log  V 
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which    linearizes   equation   (4): 


V      =  pV 


Solutions    of   the    linear   heat    equation   can   be   explicitly  expressed  by  an 
Integral  formula.      This  formula  for  V   can  be   used   to  express   u: 

(6)  u(x,  t  ;e  )    =  /  ^-^  exp    [    -  _)    dy    /  /    exp   [  -  -)    dy 
where 

(7)  D  =  U(y)    +^^^^  .  U(y)    =     /   u(z)   dz    , 
and    u(z)   is    the   initial  value.    Letting  e   ■»■  0   in   (6)   we    get 


_  _   * 

(8)  u(x,t)    =  lim  u(x,t;e)    =  ^^I^ 

G>0  ^ 


where   y  (x,t)   is   that   value   of   y  which  minimizes  D(x,y,t).   The 
function  y   is  defined  for  almost  all  x,t. 

The  key  step  in  the  method  of  Hopf  and  Cole  sketched  above  is  the 
transformation  (5)  that  linearizes  (4).  Gardner,  Green,  Kruskal  and 
Mlura  have  discovered  [3]  that  the  KdV  equation  can  be  linearized  by 
the  so-called  scattering  transformation.  In  this  study  we  use  this 
linearization  to  trace  the  dependence  of  solutions  on  e  and  to 
determine  their  limit  as  e  +  0.  The  limit  in  this  case  is,  for  t  large 
enough,  only  a  weak  limit.  The  initial  functions  we  consider  are 
assumed  to  tend  to  zero  fast  enough  as  |x|  -»•  <» ,  and  to  be  nonposltlve. 
The  latter  condition  is  used  essentially  in  the  analysis;  it  has   been 
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shown      in      the      dissertation      of      Venakides      [13]      how      to      remove   this 
restriction. 

The  organization  of  the   paper   is   as   follows: 

In  Section  1  we  solve  the  direct  scattering  problem  for  given 
initial  data,  assumed  for  simplicity  to  have  a  single  local  minimum. 
The  direct  problem  is  solved  asymptotically;  we  then  replace  the  exact 
initial  data  u(x)  by  an  approximate  one  u(x;e)  whose  exact  scattering 
data  are  equal  to  the  approximate  data  of  u(x).  As  shown  subsequently, 
lim  u(x;e)    =  u(x)    in   the  L2(R)   sense. 

In  Section  2  we  use  the  Kay-Moses  explicit  solution  of  the 
reflectionless   inverse  problem,   and  carry  out   the   limit  e   ■»■   0.   We      show 

that      lim     u(x,t;e)    =u(x,t)    exists    in  the   sense   of  weak,  convergence   in 

£->-0 
L2(R)  with   respect    to  x,  and    that    the   weak    limit  u   can  be    described   as 


(9)  u  =  a^^Q*. 


The  function  Q   (x,t)   is   determined  by  solving      a      quadratic     programing 


problem,    i.e. 


,*, 


(10)  Q   (x,t)    =     min     Q(\});x,t)    . 

Here  Q(\|;  ;x,t)  is  a  quadratic  functional  of  \|; ,  which  depends  linearly  on 
the  parameters  x,t.   The  function  ij)  is  determined  by  the  initial  data: 


(11)  4,(n)=Re/  V  ,,„  dy 

(-u(y)-n2)l/2 
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We    also   show  that    u^(x,t;e)    leads,    in  an  appropriate  weak  sense,      to     a 
limit. 

In  Section  3  we  show  that  Q  is  continuous  in  a  weak  sequential 
topology,  and  that  the  space  of  admissible  functions  is  compact  in  that 
topology.  We  further  show  that  Q  is  a  strictly  convex  function;  since 
the  admissible  functions  form  a  convex  set,  this  implies  not  only  that 
the  minimum  of  Q  is  taken  on  at  a  unique  function,  but  that  this 
function  is    the   only  one  which  satisfies    variational  conditions. 

The  variational  conditions  are  then  converted  to  a  Rieraann-Hilbert 
problem,  i.e.  to  the  problem  of  determining  an  analytic  function  of 
class  HP  in  the  upper  half  plane  whose  real  and  imaginary  parts  are 
prescribed   on  alternate    intervals    of    the   real  axis. 

In  Section  4  we  solve  the  Riemann-Hilbert  problem  for  all  times  t 
that  do  not  exceed  t,  called  the  break-time,  defined  as  the  time  beyond 
which   the   initial  value  problem 

(12)  u^   -   6  uu^  =  0    ,      u(x,0)    =  u(x) 

has  no  solution;  equation  (12)  is  obtained  by  setting  e  =  0  in  the  KdV 
equation  (1).  We  show  that  for  t  <  t^  ,  u(x,t;e)  tends  to  the  solution 
of  (10),  and  that  this  convergence  is  strong  convergence  in  L2(R)  with 
respect   to  x. 

In  Section  5  we  solve  the  Rieraann-Hilbert  problem  for  any  x  and  t 
and  verify  that  the  solution  so  obtained  satisfies  the  variational 
conditions  and  therefore  solves  the  minimum  problem.  The  solution  u 
determined  from  (9)  has  the  form  suggested  by  Whitham  [14]  and  more 
generally  by   Flaschka,   Forest   and  McLaughlin    [2]. 
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In  Section   6  we  determine   the  asymptotic  behavior     of     u(x,t)     for 
large   t,   For  0  <  x/t    <  4m^,    where 

m     =     -  mln  u(x)    , 
we  have 


u(x,t)   =    -  - —  (jjC/xMt)    ; 

ZiT  t 


here  (j)    is    the   function  defined   in  (11). 
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1.  Asymptotic  Analysis  of  the  Direct  Scattering  Problem. 

In  a  remarkable  paper  [3],  Gardner,  Greene,  Kruskal  and  Miura 
(GGKM)  have  shown  how  to  solve  the  Initial  value  problem  for  the  KdV 
equation 

(1.1)  u^  -  euu^+e\^^=   0 

for  Initial  data   u(x)    that   tends   to   zero   so   fast   as    I  xl    >  <»   that 

(1.2)  /    (1    +  x2)     |u(x)  I    dx  <  "    . 

Their  method  associates    to  each   solution  of  KdV  a    one-parameter     family 
of  Schroedinger  operators    ^(t): 

(1.3)  ^(t)    =  -  c^l  +  u(x,t). 

Recall     that    for  a  Schroedinger  operator    ^  whose   potential  u  satisfies 
(1.2)    one   can  define    scattering  data;    these   consist    of 

(I)  The    reflection   coefficient   R(k), 

(II)  The   eigenvalues   "n^,    n  =    1,...,N   in  the    point   spectrum  of 
(rij^  >    0  by   convention), 

(ill)    The   normlng   constants    c^  associated  with   the  eigenf unctions 

We  recall  the  definition  of  the  normlng  constants:  Let  f  be  the 
eigenf  unction  of  (^with   eigenvalue  -n    : 

(1.4)  ^f   =   -e2f^^  +  ^f   =  _^2f^ 
normalized  by 
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as 

(1.4)'  /  f2  dx  =  1. 

It  is  easy  to  show,  using  condition  (1.2),  that  as  |x|  -»•<»,  f(x)  decays 
exponentially.   More  precisely,  there  is  a  constant  c  such  that 

(1.5)  f(x)  =  c  e"^^/^  as  X-.-  -H». 

This    constant   c,    chosen  to  be   positive,    is   the  normlng  constant. 

The  operation  relating  the  potential  u  of  to  the  scattering  data 
{R(k.),n  ,c  }  is  called  the  scattering  transform.  It  can  be  inverted 
with  the  aid  of  the  Gelfand-Levitan-Marchenko  (GLM)  equation  [4].  GGKM 
have  shown  that  if  the  potential  u  varies  with  a  parameter  t  so  that 
the  KdV  equation  (1.1)  is  satisfied,  then  the  scattering  data  vary  in  a 
particularly  simple   manner  with   t: 

(i)        R(k,t)    =  R(k)  e^i'^  ^'^ 

(1.6)  (ii)      ri^      is    independent    of   t 
(iii)  c^(t)    =  c^   e^'^'^^'^ 

Formulas  (1.6)   have   been  derived  from   those      given     by     GGKM,      who 
treat   the  case  e   =    1,    by    the   simple   rescaling 

(1.7)  u(x,t)    =  v(x/e,    t/e). 

A  similar  rescaling  has  been  applied  at  the  end  of  this  section  to  the 
fomulas  of  Kay  and  Moses  for  the  inverse  scattering  transform. 

The  GGKM  solution  of  the  initial  value  problem  for  KdV  is  to  take 
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the     scattering     transform     of      the      initial     data,    then  apply    (1.6)    to 
determine   the    scattering  data   at   time   t,   and    then   invert   the   scattering 
transform. 

In  this  paper  we  shall  study  initial  data  which  satisfy  two 
further   conditions   in  addition   to    (1.2): 

(1.8)  u(x)    <  0, 

(1.8)'   u(x)  is  C  and  has  only  a  finite  number  of  critical  points. 

Condition  (1.2)  is  there  to  make  the  machinery  of  scattering 
theory  work  smoothly.  Formula  (1.8)  is  a  genuine  restriction  that 
allows  us  to  neglect  reflection.  Condition  (1.8)'  is  there  for 
technical  reasons.  In  fact  in  this  study  we  shall  assume  that  u(x)  has 
a  single  critical  point;  we  denote  by  Xq  its  location.  The  changes 
necessary  to  handle  the  more  general  case  (1.8)'  are  not  hard  to  make. 

We  study  now  the  asymptotic  behavior  of  the  scattering  data  as  e  > 
0.  The  distribution  of  the  eigenvalues  is  governed  by 

Theorem  1.1  (Weyl's  law):  Let 


^  =  -  eh^   +  u 


be  a   Schroedinger  operator  whose   potential   satisfies    (1.2).      Denote      by 
N(e )   the   number   of   eigenvalues   of  ^;    then 


Re 


(1.9)  N(e)  =  _/  •-u(y)     dy. 

ire 


2 
More  generally,  the  number  N(e,n)  of  eigenvalues  less  than  -n   is 


-9- 


(1.10)  N(e,n)  =  -^J   /-u(y)-n2 

ire   ' 


dy 


Suppose     u  satisfies   (1.8)   and  has   a   single  rainimum.      Then  for  -n^ 
>   min  u  there   are    two  functions    x^(y)   and  x_(n),    defined  by 

(1.11)  u(x_)    =  u(x+)    =  -  n^    ,        x_   <  x+ 

Clearly, 

u(x)    <  -  n        for     x_  <   X   <  x^.    , 

so   that    (1.10)    can  be   written  as 


(1.12) 


N(e,n)  =   —  $(n) 


where 
(1.13) 


X4.(ti) 


$(n)    =         /  /-u(y)-n^   dy    , 

x_(n) 


see  Figure   1; 


X   (n) 


We    shall   from  new   on  assume    that 
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(1.13)'  u(xq)    =  mln  u(x)   =   -    1    . 

Then  the    domain  of  the    functions    x_,   x^  and  $    is  0  _<  n  _<   !• 

We    turn  next      to     the      asymptotic      determination  of      the      norraing 

constants.        We     shall     use      a   crude  WKB   method,   i.e.  we  represent   the 
eigenf unction   f   in    (1.4)   as 

(1.14)  f(x)      =     c    e"Q(^^/^   . 

Substituting   this    into   (1.4)    gives 

"  ^x  "''  ^    ^xx  +  u  +  n^      =     0     . 
We    new    set  e    =    0  in   this   relation  and   define  9    as    the   solution   of 

which  yields 

(1.14)'     e(Ti,x)   =nx  +  /  n  -   (n2-Ki(y))^/2  dy 

X 

We  shall  show  below,  see  inequality  (1.26),  that 


/  v'-u(y)  dy  <  " .  ,  t '•• 

From  this    it    follows   easily    that   the      integral      (1.14)      converges,      and 
that   as   X  ->■  00 ,  '  "'v    ..:., 


(1.15)  e(n,x)   =  nx  +  0(1). 
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As   X  decreases  from  -H»>,  6  is  real  and  decreases  until  x  reaches  x,  (n  ) , 
We  denote  this  minimum  value  of  6  by 

oo 

(1.16)  e_^(Ti)  =Tix^(n)  +   /  n  -  (n^  +u(y))l/2  jy  . 

x_|_(n) 

We  now  define  ' 


(1.17)  c(n)  =e'+^^>/^ 


and  define  f  to  be  given  by  f omnia  (1.14),  with  9  and  c  defined  by 
(1.14)'  and  (1.17).  We  see  that  f  satisfies  the  eigenvalue  equation 
(1.4)  crudely,  and  the  normalization  condition  (1.4)'  just  as  crudely. 
It  follows  from  (1.15)  that  f  satisfies  the  asymptotic  relation  (1.5). 

If  we  apply  the  WKB  method,  crude  or  refined,  to  calculate  the 
reflection  coefficient,  we  get  R(k)  =  0. 

Definition.  We  define  the  modified  initial  function  u(x;e)  as  the 
function  whose  scattering  data  are: 

(i)    R(k)  E  0 

(ii)   r\      defined  by 
(1.18)  <I.(t1j^)  =  (n  -  1/2)  e^,    n  =  1 N(e), 

N(c)  =  [  J_$(0)] 

(iii)   c^  =  cdi^),    n=  l,...,N(e), 

where  [•]  is  the  "greatest  integer  less  than"  function,  * (n )  is  the 
function  defined  by  (1.13)  and  c(n )  is  the  function  defined  by   (1.17) 
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and      (1.16).        Definition      (1.18)   ..    of      n       is      the     WKB      result   and   is 

11      n 

consistent  with  Theorem  1.1. 

Theorem  1.2.  The  function  u(x;e),  defined  above,  approximates  the 
prescribed  initial  data  u  in  the  following  sense: 

L  -  lim  u(x;e)  =  u(x) 

Proof    of    this   theorem  will  be  presented   in  Section   4. 
Combining      (1.18)    and    (1.6)    gives    the   scattering  data    of  u(x,t;e), 
the   solution  of  KdV  with   initial   data  u(x;e): 

(4Ti^t  +  e  .(n))/e 

(1.19)  c^(t)    =e       '  +  .      n    =  n^. 

On  account  of  (1.18).  ,  the  functions  u(x,t;E)  are  ref lectionless 
potentials.  For  these,  Kay  and  Moses  have  given  a  simple  solution  of 
the  GLM  equation  which  leads  to  the  following  explicit  expression  of 
the  potential   u  in   terms   of    its    scattering  data: 

(1.20)  u(x,t;e)    =  92w(x,t;e) 
where 

(1.20)'  W(x,t;e)    =  -2e  ^  log  det(I  +  G(x,t;e)). 

Here  G  is    the  matrix 


(1.21)  -  '^ 


G   =  g{- c„c„l 


n         'm 
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Uslng   (1.19)   we   can  rewrite    this   as 

(1.21)'  G  =  ef     ^"^    ]     , 

where 

f^     ^^^  r  ^  -a(nn.X,t)/e 

(1.22)  gj^(x,t;e)    =  e  " 
Here  , 

(1.23)  a(n,x,t)   =nx  -  -^ ^t  -  e^(n)  , 

where  9^(n )  is  defined  in  (1.16). 

We  define  now  some  further  auxiliary  functions  that  will  be  of  use 
in  Section  2,  and  prove  a  few  of  their  properties: 

x_(n) 

(1.24)  e_(n)  =nx_(n)  -   /   n  -  (n2-Hi(y))  ^/^  ^^  . 

we  denote  the  derivative  of  *,  defined  in  (1.13),  by  -(j) : 

--     c    .  x^(ri) 

(1.25)  <i,(n)  =  -i-$(n)  =   /  \  ,,„  dy  . 

'^      x_(ti)  (-u(y)-n2)l/2 

Lemma  1.3.   (fi  is  a  nonnegative  function  which  belongs  to  L  (0,1). 

Proof:  This  follows  from  the  fact  that  $  is  a  decreasing  function 
of  n  which  is  bounded.  The  first  assertion  is  obvious;  to  verify  the 
second  we  state 
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(1.26) 


<   [    /     Iu(y)  1    (l+y2)   dy  /   _i_  dy]  ^^^   < 


l+y2 


Here  we   used   the  Schwarz    inequality,    and   assumption   (1.2). 

Lemma    1.4.      The  functions   9^  and  e_    ,    defined  by   (1.16)   and    (1.24) 
respectively,    are   continuous  and   satisy 

e+(0)     =    e_(o)     =    0  , 

(1.27) 

e_(n  )     <  nxQ  <  9+(n)     for  0  <  n   <  1   . 


Proof:  It  follows  directly  from  the  definitions  of  9^  and  9_  that 
they  are  continuous  for  n  ^  0.  To  show  contiuity,  and  (1.27),  at  n  =  0 
we  argue  as  follows: 

By  assumption  Iu(x)  |  is  a  decreasing  function  of  x  for  x  large 
enough.   So  for  large  x 

X 

(1.28)  2  /   (-u(y))^/2  >  ^(_^(^))l/2  ^ 

x/2 

Since  according    to    (1.26) 

oo 

/   /-u(y)   dy    <  », 
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it    follows   from   (1.28)    that 


(1.29)  lira  X  •-u(x)    =   0. 

x-*«> 


According  to  (1.11),  -u(x^)  =  n^;  from  this  and  (1.29)  we  conclude  that 


(1.29)'  lim  x+(n)n  =  0. 


This  shows  that  the  first  term  in  the  definition  (1.16)  of  e+(ri  )   tends 
to  0  as  ri  ->■  0. 

We  show  now  that  the  same  is  true  of  the  second  term.   Introducing 
\i    as  new  variable  of  integration  by 


u(y)  =  -  u^  ,   y  =  x^(u), 
we  have,  after  integrating  by  parts 

/   n  -  (n^-hi(y))^/2  dy  =  -  /  n-(n2-p2)^^^  —^  dy 
x+(n )  0  ^^ 

1      mx+(m) 
'^  J    — n — o    1  /o   dp    -  nx+(ri)    <  const.   max{  ux4(u  )}  . 
0  (n2ni2)i/2  ■      y<  n 


It    follows    from   (1.29)'   that    the  above    tends      to     0     as     n      -»■      0.      This 
completes    the  proof   of  Lemma   1.4. 
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2.    Asymptotic  analysis    of   the   Inverse   scattering   problem. 

In  this  section  we  carry  out  an  asymptotic  analysis  as  e  -►  0  of 
u(x,t;e),   given  by   formula   (1.20): 

(2.1)  u(x,t;e)    =  8^W(x,t;e), 
where  W   Is   defined  by   (1.20)': 

(2.2)  W(x,t;e)    =  -28^  log  det(I+  G(x,t;e)) 

andG(x,t;e)   Is   given  by   fomula   (1.21). 

We  expand  det(I  +  G)  by  grouping  together  all  terms  which  only 
contain  factors  from  G  with  indices  from  a  set  S.  The  resulting  formula 
is 


(2.3)  detd  +  G)   =  I    det  Gg    , 

S 


where  Gg  is  the  principal  minor  of  G  obtained  by  striking  out  all  rows 
and  columns  whose  index  lies  outside  S.  In  the  sum  (2.3),  S  ranges  over 
all  subsets  of  the  N  indices  and  we  take  det  Gg  =  1  when  S  is  the  null 
set . 

We  can  factor  G,  given  by  (1.21)',  as 


G  =  eD( — 1 — ]l>    , 


'%-*^ra 


where  D  is  the  diagonal  matrix  with  entries  gi,...,gfj.  Since  the 
principal  minors  Gg  of  G  have  the  same  form  as  G,  a  similar  expression 
holds  for  them.  The  matrix  in  the  middle  is  of  Cauchy  type;  its 
determinant  is  given  by 
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n       in^-^^l   /    n     (n^+n^)  , 

Sx  '  S  Sx  S 


where  Sx 'S  means  SxS  minus   the      diagonal.         Since      det     D     =     II      g        we 

S 

conclude   that 


detGs  =  [^)isi  nfn   n      1^1 

S      "   Sx'S        "     ™ 


Using    (1.22)    we   can  write    this    as   follows 


(2.4)  det  Gg  =  (L)ISI    exp   {-q^/2eh    , 


where 


(2.5)      Qs(x,t;e)    =  %   I    a(nn,x,t)    -  e^     I        iog(!l^^)  2  +  2  e  2  j;    i^g  ^^. 

S  Sx'S  '""^"^  S 


Using    (2.4)    to  express   (2.3)    gives 


(2.6)  det   (I  +  G)    =  I        [i)     IS"    exp   (-Qs/2e2)    . 

all  S 


We   shall   show   that   for  small  e   this  sum  is  dominated  by  its 
largest  term.   Define 

(2.7)  Q*(x,t;e)  =  min  {  Q3(x,  t  ;e  )}  . 

S 

Since  all   terras    in  (2.6)   are   positive,    and   since      there     are     2^^^^  ^      of 
them,    we    get    the   following    two  sided    inequality 
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(^]N  ^-Q*/2e2  <  ^^^(j  +  G)  i  2N  e-Q*/2E 


Taking  the   logarithm  and   multiplying  by   2z      gives 


2tH  log   d.)    -  Q*   <   2e2  log  det   (I  +  G)    <  2£  ^N  log  2    -  Q* 


The  center    term  is   just  -W     as      defined     by      (2.2)      so     we      obtain     the 
estimate 


|W(x,t;e)   -  Q*(x,t;£)|    <  Iz'^nU)   log   (i). 

—  e 


Since  according  to  (1.9),  N(e )  =  0(e   ),  we  have  proved 
Theorem  2.1. 


(2.8)       lim  [W(x,t;e)  -  Q*(x,t;e)]  =  0, 
e-»-0 

uniformly  in  x  and  t. 


We  shall  show  that  the  limit  of  Q  (x,t;e)  exists  and  can,  just  as 
each  Q  (x,t;e),  be  characterized  as  a  solution  of  a  minimum  problem. 
To  find  this  rainlraum  problem  we  rewrite  Qg(x,t;e)  given  by  (2.5)  as  a 
Stieltjes    Integral. 

For  each    S  we   define  i|)g(n,e)  as    the   distribution 


(2.9)  «^s(ri,e)    =  e   ^    I    ^{r\-^^) 

S 

and  rewrite  (2.5)   as 
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.    1 
(2.10)     Qs(x,t;e)  =_/  a(n,x,t)  iPgCn.e)  dn 

0 

1    ^    ^ 
-  —  /  /   log  (in^)2  ^g(p  ,£  )x  >g(n  ,e  )  dp   dn 

TT^  0   0  '^"^ 

2e    ^ 
+  —  /   log  n  'i'cCn.e)  dn   ; 

""    0 


here  \(/g(u  ,e  )  x'  \|>g(n,£)  is   the  product   of   the  distributions     minus      the 
diagonal  terms, 

(2.11)  ^^iv,e)  ^'  ^s^"^'^^  '^^^    I      «(p  -  n^,)  6(n  -n^). 

Sx'S 
We   define  ()>  (n  ,e  )   by 

N(e) 

(2.12)  <|)(n,e)    =£71        j;      6(n   -  nn)    ; 

n=l 

Clearly 

(2.13)  -.     ,  0  <i|;s(n,e)   dn    <(t)(n,e)   dn 


as  measures,    for  every   S. 

Since   the  n      satisfy  Weyl's    law  (Theorem  1.1),    it   follows   that 

(2.14)   lim  (|,(n,e)  dn  =  <|)  (n  )  dn, 
e+0 

,  '  -     lim  (fi  (y  ,e  )  X  '  (!>  (n  ,e  )  dpdn  =  (|)  (vi  )  (J)  (n  )  dy  dn  , 
e+0 

in  the  sense  of  weak  sequential  convergence  of  measures,  where  4"  (ti  )  is 
defined  by  (1.25).  These  considerations  suggest  a  minimum  problem 
satisfied  by  the  limit  Q*(x,t;e): 
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Theorem  2.2.  . 

(2.15)  lim  Q*(x,t;e)    =  Q*(x,t) 

E+0 

uniformly    on   compact   subsets   of   x,t,    where 

(2.16)  Q*(x,t)    =  rain  {Q(i});x,t):   <|;  e    A}, 

the  admissible  set  A  being  all  Lebesgue   measurable   functions  \|;    on   [0,1] 
satisfying 

(2.17)  0   <  4;(n)  _<  ([.(n)     ,         _ 
and  Q(ij;;x,t)    denotes    the  quadratic   form 


4 ; 

(2.18)  Q(i|;;x,t)    =-/    a(n,x,t)   i|;(ri)    dn        , 

IT 

0 

11 

-  _L/   r  log  (^)2  ,|,(u)  ,|;(n)  dy   dn 


The  proof  is  based  on  a  series  of  lemmas  and  theorems. 

Let  A  denote  the  set  of  distributions  i|;g(Ti,e)   defined  by   (2.9) 
for  all  possible  sets  of  indices  s  and  e  >  0  , 


(2.19)         A^  =  {4is(n,e):  S,  e}  . 

To  the  elements  ij;  of  both  the  admissible  sets  A^  and  A  we  assign 
positive  measures  \J;  dn  on  [0,1],  If  ij^g  e  A^  then  from  (2.13),  (2.12), 
and  (1.18)^^ 

.  ,■■..  1   .  -.■       .  p   ^-  -  ■ 


(2.20)       /  ^s(Ti,e)  dn  ^  /  (t)(n,e)  dn  =  e  TT  N(e)  <j   /-u(y)  dy  , 
0  0  -» 


while  if  ij;  e  A  then  from  (2.17)  and  (1.25) 
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(2.21)         /  ,Kn)  dn  <  /  4>(n)  dn  =  /  /-u(y)  dy 

0  0  -co 


Thus  we  conclude  from  Lemma  1.3: 

Lemma  2.3.  The  total  variations  of  the  measures  \|)  dp  ,  i|;  in  A  and 
A  are  uniformly  bounded. 

Definition.  A  sequence  of  distributions  ij;.  is  said  to  w-converge 
to  i|;  if  the  corresponding  measures  ij/.  (n)  dn  converge  to  ij)  (n  )dri  in  the 
sense  of  weak  sequential  convergence  of  measures,  that  is  if 


(2.22)  lira  (x,tk)  =  (X.'l') 

k-H» 


for   every    continuous   function  x*      Here    the   parentheses   denote   the  usual 
duality   on   [0,1] 

1 
(2.23)        (x,^\>)   =  /  x(n)  *(n)  dn  . 

0 

We  shall  abbreviate  (2.22)  as  w-convergence  and  write 


w-lim  ij;,  =  \|; . 
k-xjo 


The  next  lemma  characterizes  the  admissible  set  A  as  w-limits  of 
the   distributions    in  A    . 

Lemma  2.4.  (a)  Let  \j;,  denote  a  sequence  of  elements  in  A^  ,  i|'{^(n ) 
-  'I's    (n,ej^),    such   that   as  Ej^  -»•   0, 


(2.24)  w-lim  ^^(t))    =  ^(n); 

k-x» 


then  \|;  belongs  to  A. 

(b)  Conversely,  for  each  i|j  in  A  we  can  choose  S(e)  so  that 
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(2.25)  w-lira     'l'S(e  )('l '^  >    =  *(^)' 

(c) 


(2.26)  w-lira     ((>  (n  ,e  )    =  <J)  (n  )    . 

£-.-0 


Proof:    By  (2.13)    the  ^l)^  satisfy 

0  _<  ^]^(r] )   dT]    <  <^(T],e^^)   dn    . 

Since  inequalities  between  measures  are  preserved  for  their  weak 
limits,  it  follows  from  (2.24),  and  (2.26)  that 

0  ^  ij;  (n  )  dn  <  (()  (n )  dn . 

From  this  we  conclude  by  the  Radon-Nikodym  theorem  that  ij;  is  a 
measurable   function  and    that    (2.17)    is    satisfied;    thus  \(;  e    A. 

Part(b)  is    straightforward  and  part   (c)    is    contaied   in   (2.14). 

The  main  step  in  proving  Theorem  2.2  is  in  relating  Qc(x,t;G)  as 
given   in  (2.10)    to  Q((|;;x,t)   as    given   by   (2.18). 

Theorem  2.5.  Denote  by  \])y^  a  sequence  of  elements  in  A^.  ,  i^]^(n  )  = 
tg    (n,e^^),    such   that   as  e^^  ^-   0, 


w-lira  ^\)y.ir\)    =  4)(n)      ; 
k-x=° 


then 


(2.27)  lira  Qg    (x,t;ei,)    =  Q(,p;x,t) 

k->0        '^ 
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for   every  x  and   t. 

The  mild  technical  difficulties  in  the  proof  arise  in  handling  the 
second  and  third  terras  of  the  right  hand  side  of  (2.10).  This  will  be 
accomplished  with  the  aid  of  lemmas,  and  the  introduction  of  the 
operator  L: 


(2.28)        L4;(n)    =  J-    /    log   (115^)2  i(,(u)   du 


2ti     ■'q  n+w 


Lemma  2.6.  Lji  is  a  continuous  function  on  [0,")  that  tends  to 
zero  as  n   ^  °°  •      In    the   interval    [0,1] 

(2.29)  -L(j,(n)  =  e+(n)  -  9_(n)  . 

Proof;  Since  it  was  shown  in  Lemma  1.4  tha  9+ (n )  are  continuous 
functions,  the  continuity  of  LcJ)  in  [0,1]  follows  from  (2.29).  The 
continuity  of  L<J)  on  [I,")  and  behavior  at  infinity  can  be  seen  directly 
from  formila  (2.29). 

Foritula  (2.29)  will  be  derived  at  the  end  of  Section  3. 

Corollary  2.7. 

(a)  For  every  \p    in  the  admissible  set  A,  L\j;  is  continuous  and 

(2.29)'     L(t)  (n  )  _<  Li|;  (n  )  _<  0  for  n  >  0. 

(b)  L\l)   is  an  odd  function  of  n  • 

Proof:  (b)  is  obvious  from  (2.28)  so  we  restrict  our  attention  to 
n  >  0.   Inequality   (2.29)  '  of   (a)   follows   from  the   fact   that  by 
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ineqtiality   (2.17)   admissible   functions  are  nonnegative,  and  that  the 
kernel  of  the  integral  operator  L  is  negative. 

To  show  the  continuity  of  Li|;  at  n  ,  n  >^  0  we  take  any  sequence  {lu} 
converging  to  n  with  n^  >^  0.  The  integrands  of  Li|;(n{^)  appearing  in 
(2.28)  converges  almost  everywhere  to  the  integrand  of  Li);  (n  )  and  are 
bounded  below  by  the  integrands  of  L<|((n^).  Since,  by  the  continuity  of 
L(j)  ,  we  know  Lij)(n.  )  converges  to  L(t)(ri),  we  conclude  by  a  dominated 
convergence  argument  that  L\|)(ri|^)  converges  to  Lij;(n).  The  continuity  of 
Lij;  follows,  completing  the  proof. 

We  use  the  operator  L  defined  by  (2.28)  in  conjunction  with  the 
duality  notation  of  (2.23)  to  rewrite  the  functional  Q(^;x,t)  of  (2.18) 
as 

(2.30)  Q(i^;x,t)  =  1  [(a(n,x,t),t|;)  -  1.  (L,|;  ,,|; )]  ; 

note  that  by  Lemma  1.4,  a(ri,x,t)  given  by  (1.23)  is  continuous.  In 
particular  this  shows  Q(i);,x,t)  is  finite  for  all  \j;  in  the  admissible 
set  A. 

Preparatory  to  proving  Theorem  2,5,  we  break  up  the  functionals 
Q(ij;)  and  Qg  into  two  parts  depending  on  a  parameter  B,  which  eventually 
tends  to  infinity: 

(2.31)  Q(,|;)  =  qS(^)  +  Qg(^)  ^ 
and 

(2.32)  Qs   =Qi+QB,g  ... 
Here,  setting 
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'<"-'    ■  -ii"«(S''  • 


and 

(2.33)  £^(n,w)      =     mln  {£(n,y),B}     , 

we    define 


,1  ,11 

(2.34)     qB(;j,)  =lj  a(n)'|^(n)dn  +4-11  £^(n,p)<|'(n)tKy)  dn   dy 


,1  11 

(2.35)     q|    =    -/  a(n)  'l's(n)  dn  +  -  J  /  £^(n,y)  ^-s^i)  'l's(w)  *i   ^^ 
'^    0  ^00 


Lemma  2.8.      (a)    If    the   sequence  \|;^  =  ifig    (n  ,G]^)   satisfies 


w-lim  \|;j^(ri)    =  t    » 
then 


^'™Qs  (^k>  =  Q^(*>  • 

k-HX.  ^k  ^ 


(b)   For   all  ij;    in  A 

(2.36)  \QQii>)\    <  6(B)      , 
where 

(2.37)  11m  6(B)    =   0 

B+0 


(c)   For   all  ii{e)   in  A^ 


i^s.bI    <  ^^^^  ^  ^(^^ 
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IQs,bI    <  '5(B)   +  Y(e)      , 
where 

lim  y  (e  )    =  0  . 
e+0 

Proof:  Since  both  functions  a(n )  and  £"(n,u)  are  continuous,  part 
(a)  is  a  direct  consquence  of  the  definition  of  weak  convergence 
(2.22). 

To    prove   (2.31)    and    (2.34)  we   write 


2 
%(.^)   =-//      (£  (n  ,P  )-B )  4- (n  )  If  (u  )  *!   d;i 
"^    £>B 


Since  for  admissible  \\) , 

it  follows  that 

IQbC'I')!  <  QbC*)   • 

This   proves      (2.36)     with     6(B)      =     Qb((},  )    ;      (2.37)      follows      from     the 
observation   that     Q(()))    <  °°    . 

For   the   proof   of   (c)  we    refer   to    [9    ];   we    remark   that    it    relies   on 
the   choice   (1.18)^.      of   the   n     to    show  separately        that 


Jyl  ^('^n.nn.)    -B    <   6(B)    . 


and   that 

N 
2 


el    log  '^  n     <  y(e) 
1 
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In  showing  that  Y(e)->-  Oase  ■*■   0  one  has  to  use  the  fact   that   for  a 
potential  u  that  satisfies  (1.2), 

(2.38)  /  (-u(x))l/2  log  (-u(x))  dx  <». 

Theorem  2.5  is  an  immediate  consequence  of  Lemma  2.8. 

The  pointwise  limit  (2.69)  of  Theorem  2.5  is  uniform  on  compact  subsets 
of  x,t  space.   This  follows  from 

Lemma  2.9.   The  family  {Qg(x,t;e):  S,e}  of  linear  functions  of 
X  and  t  defined  in  (2.10)  is  equicontinuous  and  equibounded  on  compact 
subsets  of  X, t  space. 

Proof:   Since  the  functions  Qg(x,t;e)  are  linear,  it  suffices  to  show 
the  coefficients  of  x,  t  and  the  quantity  Qg(0,0;e)  are  uniformly 
bounded.   The  coefficients  of  x  and  t  are  uniformly  bounded  since 


S 

e  I  Hn  <  e  N   , 
S 


respectively  and  by  (1.9)  eN  is  bounded.  A  uniform  bound  for  Qg(0,0;e) 
will  be  demonstrated  by  considering  the  three  terms  of  the  right  hand 
side  of  (2.10)  separately.  According  to  Lemma  1.4,  e^(n)  is  continuous 
and  therefore  bounded,  say  6^(0)  _<M.  The  first  term  of  Qs(0,0;e)  is 
then  uniformly    bounded   since 


e   I    e+(nn)  i  e   N   M    , 
S 
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and  by    (1.9)   eN  is    bounded.      The  uniform  bounds  of    the   second  and    third 
terras    of  Qg(0,0;e)   follow  directly    from  Lemma  2.8. 

We   now  turn   to  Theorem  2.2. 

Proof:    Formula    (2.15)    of   Theorem  2.2  will   follow     from     two      facts 
proved  below: 


(2.39)  11m   sup     Q*(x,t;e)  _<  inf   {Q(il;;x,t):   \|;  e    A}    , 


e   ■►    0 
and 


(2.40)  lira   inf  Q*(x,t;e)    =  Q(il;*;x,t)      , 


e   -^   0 


for   some  4;      in   the   admissible   set  A.    Comparing    (2.39)   and    (2.40)  we    see 
that 


(2.41)  Q(i|;*;x,t)    =  inf  {Q(4);x,t):   i>  e    k] 


Thus  the  minimum  asserted   in  the   definition   (2.16)   of  Q  (x,t)   is 
attained, 


(2.42)        Q*(x,t)   =  Q(^*;x,t)  , 


and  (2.15)  follows. 

To      prove   (2.39),    let  ip    be  any   element   of    the  admissible   set  A.    By 
Lemma  2.4,    there   exists  S(g )   such   that 


(2.43)  w-lim     iPs(e)<n,e)    =  t  (n )      . 

e+0 


l^t   {e^}    be   a   sequence   such   that  e^^  ->•   0  and 
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(2.44)        lira  Q*(x,t;e^^)    =  11m  sup  Q*(x,  t;e  ) 
k.-H»  e   +   0 


By  definition   (2.7)   and  Lemma  2.9,    theQ*(x,t;e)   are      the     minima      over 
subfamilies   of   an  equlbounded   family    of    functions   and    so  are    themselves 
equibounded ,    thus    insuring    the    limits    of    (2.44)   are   finite. 
By    (2.7) 


Q*(x,t;e^^)     _<  Qg    (x.f.E^^)    . 

K. 


Using  (2.43)    and   applying  Theorem  2.5  we    see   that    the   right   side      above 
tends    to  Q(\|;;x,t).      Using    (2.44)    on   the   left  we    obtain   in  the    limit 


(2.45)  lim  sup     Q*(x,  t  ;e  )  _<  Q(i|/ ;x,  t) 

e   ->•   0 


for   all  ip  e  A .  Since  the  left  hand  side  of  (2.45)  is  independent  of  \l) , 
we  take  the  infimura  of  the  right,  obtaining  (2.39). 

To  prove  (2.40),  let  S(e)  be  such  that  the  minimum   of   (2.7)   is 
attained, 


(2.46)         Qs(e)(x,t;e)  =  Q*(x,t;e)   . 


Let  {e^^}  be  a  sequence  such  that  e^^  ->■  0  and 


(2.47)         lim  Q*(x,t;e^)  =  lim  inf  Q*(x, t ;e ). 
k-H"  e  -•■  0 


Since  by  Lemma  2.3  the  total  variations  of  the  measures  <i's(e  )(n,e)  dn 
are  uniformly  bounded,  we  apply  the  Helly  selection  theorem  and  extract 
a  w-convergent   subsequence.    Passing   to  this  subsequence  we  know  by 
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Leraina    2.4   that  the   w-llmit  i};      lies   in   the  admissible      set      A.      Applying 
Theorem  2.5    to   this   sequence   and   using    (2.46)    we    obtain 


(2,48)  Q(^*;x,t)    =  lira  Q*(x,t;e^), 


which,    along  with  (2.47),    proves    (2.40). 

The  limit  in  (2.15)  is  uniform  on  compact  subsets  of  x  and  t  since 
by  definition  (2.7)  and  Lemma  2.9,  the  Q  (x,t;e)  are  minima  over 
subfamilies  of  an  equicontinuous  family  of  functions  and  so  are 
themselves  equicontinuous.  The  result  then  follows  from  the 
Arzela-Ascoli    theorem  and    the    proof   of    Theorem  2.2   is   complete. 

Combining  the  limit  (2.15)  of  Theorem  2.2  with  (2.8)  of  Theorem 
2.1  we    conclude    that 


(2.49)  lim  W(x,t;G)    =  Q*(x,t) 

e+O 


uniformly    on   compact   subsets   of   x  and   t. 

Let  us  denote  convergence  in  the  sense  of  distributions  of 
functions  of  x  and  t  by  d-lim.  Since  derivatives  of  a  uniformly 
convergent  sequence  of  functions  converge  in  the  distribution  sense  we 
conclude    from    (2.1)    and    (2.49) 

Theorem  2. 10.  Let  u(x,t;e)  be  the  solution  of  the  KdV  equation 
(1.5)  with    initial  data    u(x;e )    given  by   (1.18);    then 


(2.50)  d-lim  u(x,t;e)    =  u(x,t) 

e^-0 


exists,    and 
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(2.51)  u(x,t)    =  3^^  Q*(x.t) 


where  Q    (x,t)    is  defined   by    (2.16). 

Remark.  Actually  we  have  proved  more;  we  have  shown  that  u(x,t;e) 
tends  to  u(x,t)  in  the  x-dlstributlon  sense.  I.e.,  for  any  C^  function 
w   of    X, 


(2.52)  lim    (u(t;e),w)    =   (u(t),w)    ; 


furthermore  the  limit  is  uniform  over  compact  subsets  of  t.  We  will 
further   strengthen   this    result    in  Theorem  2.14. 

Exploiting  the  concept  of  convergence  in  the  distribution  sense  we 
shall   now  deduce  as   a   corollary   of    Theorem  2.10   the  following   results: 

Theorem  2. 1 1.  Let  u(x,t;e)  be  as  in  Theorem  2.10,  then  the 
d- limit 


(2.53)  d-llm  u2(x,t;e)    =  u2(x,t) 

e-»-0 


exists,    and 


(2.54)  u2(x,t)    =i.  3„tQ*(x,t). 


7  °xt^ 


Similarly 


(2.55)  d-lim   [u3(x,t;e)   +  3- e  ^u2(x,t  ;e  )  ]    =u^(x,t) 

e>0 


exists  ,    and 
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(2.56)  u^(x,t)    =  JL  9j.j.Q*(x,t) 


Proof:   We   rewrite   the  KdV  equation    (1.5)    in  conservation   form 


u,   -  (3uh^  +  t\^^=   0. 


Substituting    (2.1)    into    the    first   term  and    integrating  with   respect      to 
X  we    see    that    u(x,t;e)   satisfies 


9xtW-    3u2   +e^u^^  =   0   . 


Here  we  ascertain  that  the  constant  of  integration  is  zero  by  letting  x 
^  ».   Solving  for  u  we  obtain 


2 
(2.57)  u2  =^a^tW  +  V^xx  • 


The  first   term  on   the   right  hand   side      of      (2.57)      has     by      (2. 49)      the 

d-limit     i  8^^Q*  while  by    (2.50)    the   d-limit   of    the    last    term   in   (2.57) 
3      xu 

is  zero.      This  shows   that    (2.53)   and    (2.54)    hold. 

To  prove    (2.55)    and    (2.56)    we  multiply   the   KdV  equation      (1.5)      by 
2u   and  write   the   result    in   conservation    form: 

(u2)^   -   (4u3)^  +  e2[(u2)^^  -    3u2)^  =   0   . 

Substituting    (2.57)    into  the    first   term  and    integrating  with   respect    to 
X  we    see   that    u(x,t;e)   satisfies 
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2 


T  9ttW  +  V  "xt    -    ^-^   +  ^'(-^)xx  -    3eM   =    0 


Here,    as    before,   we  ascertain   that    the   constant    of   integration   is      zero 

o  3       9    ? 

by    letting  x  -»■  « .      Solving  for  u""  +  _  e   u^  we  obtain 

4  X 


(2.58)      u3+le2u2   =J_a^^W-H£iu,,  +^(u2) 


XX 


The      first      term     on     the      right     hand    side   of    (2.58)   has   by    (2.49)    the 

1  * 

d-limit  __  9      Q    ,    while    the    d-limits   of    the    last    two      terms      in      (2.58) 

are   by    (2.50)   and    (2.53)    zero.      This  shows   that    (2.55)   and    (2.56)    hold, 

finishing  the    proof. 

Theorem  2.10   can  be    strengthened: 

Theorem  2.12.  Let  u(x,t;e)  be  as  in  Theorem  2.10.  As  functions 
of  X,  u(»  ,t;e)  converges  to  u(«  ,t)  weakly  in  L  (R)  as  e  -»■  0  uniformly 
over   compact   subsets   of    t. 

The    proof    requires 

Lemma  2.13. 


(2.59)      lim  /    u2(x,t;e)   dx  =  /   u2(x)    dj 
e  +  O-oo  -«> 


uniformly  in  t. 

Remark.  It  is  not  surprising  that  (2.59)  is  independent  of  t, 
since  J  u  dx  is  conserved  for  solutions  of  KdV. 

Proof:  Making  use  of  the  explicit  form  of  W(x,t;e),  (2.2-2.5)  and 
the  fact  that  u^  vanishes  as  |x|  ->■«',  we  integrate  (2.57): 
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OD 

(2.60)      J   u^(x,t;e)   dx  =  i.  3^W| 


16       lall  s(^S  "^n^    ^^^  ^S  " 


e 


~^  ^all  S    "^et  Gg 


N 
16  ^  3 

-3-^      A        In    • 
n=l 


3 

Here  we   have    taken   limits    of      a     weighted      mean      of     Eg     n^     where      the 

weights,  det  Gg  ^  are  decreasing  exponential  functions  of  x,  given  by 
(2.4).  At  the  upper  limit  the  dominant  term  occurs  when  S  is  the  null 
set  while  at  the  lower  limit  the  dominant  term  occurs  when  S  consists 
of   all   indices . 

It    follows   from    (2.26)    of   Lemma    2.4   and  the    definition      (2.22)      of 
w^convergence   that 

N 
(2.61)  lim  zi,      I       nl  =    (n\<t>)    • 

e>0  n=l 

Using  the  definition  (1.25)  of  ({)  (n ) ,  interchanging  the  order  of 
integration  and  performing  the  n  integral  gives 


1  '^+(^)        4 

(n^.cf)  =/   /   ""  dy  dn 

0  x_(n)  (-u(y)-n2)l/2 


.4 


«>  /-u(y)  —   dn  dy 

(2.62)  =  /    /     (-u(y)-n2)l/2 

-00    0 


=  ^  /  u2(y)  dy 
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Putting    (2.60-2.62)    together  we   obtain    (2.59)   and   prove   the    lemma. 

Proof   of   Theorem  2.12:    By   Lemma    2.13  we    can   find   an     M     >      0      such 
that    for  all   t  and  e 


°°  1/2 

(2.63)        llu(.,t;e)ll         =(      /   u^(x,t;G)   dx)  <M    . 


Given  any  v  e    l2(R)  and  6    >    0,    we    can  find  a  C^(R)   function  w   such   that 


(2.64)  llv  -  wll    „   <        5    . 

l2     - 


We   wr  1 1  e 
(2.65)  (u(t;e),v)    =   (u(t;e),w)   +   (u(t;e),v-w)    . 

According   to   (2.52),   the  first  term  on  the  right  converges,  uniformly 
for  bounded  t,  to 

(u(t),w)  . 

The  second    term   is   by   the   Schwarz    inequality   and    the      estimates  (2.63) 

and      (2.69)      less     than  M6 .      It    follows    from  this    that    the   left  side   of 

(2.65)    tends    to  a  limit   for  every    v.    This  proves    that      u(x,t;e)  has      a 
weak      limit      as  e   -»•    0,    uniformly    for   bounded   t.    In    view   of   (2.52),    this 

weak   limit    can  be   identified  with  u.      This    shows   that    u(x,t)    is  in     L2 
with   respect    to  x,   and    proves   Theorem  2.12. 

We   end    this   section   by   remarking    that    one    can   show 


(2.66)      lim     /    u(x,t;e)    dx  =     /    u(x)    dx 
e->-0  -00  -00 
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uniformly      in      t  using  arguments   similar   to   those   that    led    to    (2.59)    of 
Lemma  2.13.       Since  u(x,t;e)  _<  0,    this   fact    can   be   used    to      replace      the 
weak     !/■   convergence    in  Theorem  2.12   by   weak   l'^  convergence    for    1   <    p   <^ 
2    just  by  using    the   obvious  LP   analogs   of    the    steps    in   the    proof. 
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3.    Analysis    of   the  Minimum  Problem. 

Theorem  2.10  gives  an  explicit  formula,  (2.51),  for  u(x,t),  the 
limit  in  the  distribution  sense  of  solutions  to  the  KdV  equation  as  the 
coefficient  of  dispersion  tends  to  zero.  The  formula  involves  the 
function  Q  (x,t)  characterized  by  a  minimum  problem,  (2.16),  where  the 
quantity  Q(\|;;x,t)  to  be  minimized,  (2.30),  is  a  quadratic  function  of  ip 
and  the  functions  v|;  of  the  admissible  set  A  are  subject  to  two  linear 
inequalities,  (2.17),  a  so-called  quadratic  programming  problem.  In 
this  section  we  begin  to  use  this  formula  to  study  a  variety  of 
properties  of  the  limit  u(x,t)  and  the  related  limits  u2(x,t)  and 
u-^(x,t)    given  by  formulas    (2.54)   and    (2.56)   of   Theorem  2.11. 

We  start    with   some   easy    observations: 

Theorem  3.1.      As  a   function  of    x  and    t,    Q     is 

(a)  <  0 

(b)  continuous 

(c)  concave 

(d)  increasing  In  x 

(e)  decreasing  in  t. 

Proof:  Since  0  e  A ,  it  follows  from  (2.16)  that  Q*(x,t)  <  Q(0;x,t) 
=  0;  this  proves  (a).  Using  the  definition  (1.23)  of  a(n,x,t)  in  the 
formula   (2.30)   for  Q(4i;x,t)   one   sees 


(3.1)        Q(i|;;x,t)    =  -  (n,«|;)x  -  —  (n^>l')t   -  -(9  +  ,;|;)    -  -   (L<|;  ,i|, ) , 

ir  IT  IT  IT 


Since   the   admissible  ^    satisfy  0  _<  i|;  ^  (J)  and  (^    belongs  to  LMO.l],  it 
follows  easily  from  (3.1)  that  {Q(ij;;x,t):  i|;  e  A}  is   an   equicontinuous 
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family      of      functions    of    x  and    t.    It    follows   that   Q*,    the    Infimum  of   an 
equicontinuous  family   of    functions,    Is   itself   continuous. 

Each  function  Q(4);x,t),  being  linear  in  x,  t,  has  properties  (c) , 
(d)   and    (e ) .      It    follows   that    so  does   Q    ,    their  infimum. 

We  draw  now  some  conclusions  from  Theorem  3.1.  From  the  concavity 
of  Q  we  deduce  that  the  matrix  of  second  derivatives  of  Q  is 
negative,  in  the  distribution  sense.  Using  expressions  (2.51),  (2.54) 
and    (2.56)    for   these  second  derivatives   yields 

Corollary   3.2. 

(3.2)  5^  +   6Ct^  +    12  T^P"  <   0, 
for  all   real  C,   x.      In    particular 

(3.3)  u   <^  0     and     "^  <   0. 

We  remark  that  there  is  no  maximum  principle  like  (3.3)  for 
solutions  of  the  KdV  equation.  That  is,  if  the  initial  values  of  a 
solution  u(x,t)  of  the  KdV  equation  satisfy  u(x,0)  <_  0,  one  cannot 
conclude  u(x,t)  <^  0.  This  may  be  seen  by  noting  that  if  u(xqjO)  =  0, 
the  value  of  the  third  x-derivative  of  u  at  Xq  could  well  be  negative. 
It  follows  then  from  the  KdV  equation  that  U|.(xo,0)  >  0,  and  so  u(xo,t) 
>  0  for  small  t  >  0.  Assumption  (1.8)  is  that  u(x)  _<  0,  thus  once  we 
show  u(x,0)  =  u(x)  then  (3.3)  can  be  interpreted  as  a  kind  of  maximum 
principle   for   the    zero   dispersion    limit   of   solutions    to  KdV. 

To  derive  further  information,  some  properties  of  the  integral 
operator  L,  defined  by  (2.28)  are  needed.  The  first  is  a  strengthening 
of  Corollary  2.7  which  asserted  that  for  every  ij;  in  the  admissible  set 
A,   Lij)  (n  )   is    a   continuous    function  on  R   that    vanishes   at    infinity. 
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Theorem  3.3.  The  set  of  functions  {L\|):  <()  e  A}  Is  an 
equlcontinuous  subset  of  Cq(r). 

Proof;  We  break,  up  the  operator  L  into  two  parts  depending  on  a 
parameter  B  which  will  tend  to  infinity: 

(3.4)  L  =  L^  +  Lg  , 

where 

1 

(3.5)  ....   lB  ^(n)  =  -/  £B(n,ii)  .|^(p)  dp   , 

0 

and  z'^   is  defied  by  (2.33).        -. 

As   B   tends   to  infinity,  lS(J)  (n  )  will  tend  to  Lcjj  (n  )  ,  monotonically 
for  all  n  e  R .  Since  these  functions  are   in  C^(R),   it   follows   from 
Dini's   theorem   that   the  convergence  is  uniform  in  n •   That  is,  if  we 
define 

(3.6)  6  (B)  =  max  |Lg<(,  |  , 
then 

(3.7)  lim  6(B)   =   0  . 

B+0 

Let  \^i    be  any  admissible  function;  since  admissibility  means  that 

. , '  .  0  ■<  iKn  )  ^  4,  (n )  , 

and   since    the  kernel   of  Lg   is    of    one   sign,    it    follows   by   (3.6)    that    for 
any  n  , 

(3.8)  |Lg,i;l    <    |Lb<J)  I  _<  6(B)    . 
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The  kernel  of  L^,  defined  in  (2.33),  is  uniformly  continuous.  It 
follows  that  L^  maps  any  set  bounded  in  the  L'^  norm  into  an 
equicontinuous  set.  Since  the  admissible  functions  are  Ll-bounded,  it 
follows  that  they  are  mapped  by  L^  into  an  equicontinuous  set.  Since 
by  (3.8),  L^p  differs  f rom  Li|)  at  most  by  6(B),  and  since  6(B)  ->■  0  as  B 
■»■  »,    Theorem  3.3   follows. 

■Recalling  the  definition  of  w-convergence  (2.22)  for  sequences  in 
A,  we  use  the  previous  result  to  prove  the  continuous  dependence  of  L 
and  Q  on  \() : 

Theorem  3 .4.      If  \()^   is   a  w-convergent   sequence   of   elements   in  A: 

w-lim  ij^u  =  ij) , 
then 
(3.9)  lim  Li|;,^(n)    =  Lil;(ri)      , 

uniformly    on  R,    and 


(3.10)  lira  Q(,|;^;x,t)    =Q((|;;x,t), 

k-x» 


uniformly  on  compact  subsets  of  x  and  t. 

Proof:  By  Theorem  3.3,  the  functions  L\|;i^  are  equicontinuous,  so  by 
the  Arzela-Ascoli  theorem  it  is  enough  to  show  the  convergence  of  (3.9) 
pointwise.   Since  the  kernel  of  L   is  continuous,  it  follows  from  the 
definition  of  w-convergence  that  for  each  n  , 
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lim  L  ij^j^(n)  =  L  iKn) 

Since  by  Theorem  3.3   L  iji   differs  at  most  by  6(B)   from  Lii,     ,    (3.9) 
follows. 

Since  the  family  of  Q(\j/;x,t)  is  equicontinuous  in  x  and  t,  it 
suffices  to  shew  the  pointwise  convergence  of  (3.10).  We  consider 
separately  the  two  terms  on  the  right  side  of  formula  (2.30)  for 
Q(\|;;x,t).  According  to  Lemma  1.4,  a(ri,x,t)  given  by  (1.22)  is  a 
continuous  function  of  n  and  so  by  the  definition  (2.22)  of 
w-convergence  .  .   , 


(3.11)  lit"  (a.'I'k)  =  (^''t') 

k-x" 


for  every  x  and  t. 

Next  we  consider  the  identity 

(3.12)   (L'^k>'l'k)  "  (.^i>,i>)      =   (L'I'k"  '^'1''  *k)  +  (.L^^),^^-^). 

The  first  term  on  the  right  is  bounded  in  absolute  value  by 

Fornula  (3.9)  shows  that  the  first  factor  tends  to  zero;  since  ^  is 
admissible  the  second  factor  is  bounded  by  ll(J)ll  ..  Thus  the  first  term 
on  the  right  of  (3.12)  vanishes  as  k  •>•  <» .  Corollary  2.7  asserts  Lij;  is 
continuous  so  the  second  term  on  the  right  of  (3.12)  tends  to  zero  as 
k  ■»■  00  by  the  w-convergence  of  ij;^  to  ij).   This  proves  that 
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(3.13)  11m  (bl'k.'I'k)    =   (Lip,*) 

k+oo 


Combining   (3.11)    and    (3.13)   according   to   formula   (2.30)    for  Q(i);;x,t)  we 
obtain  (3.10)   as   asserted   in  the   theorem. 

An    immediate   consequence   of    (3.10)    In  Theorem  3.4   Is 

Corollary  3.5.  The  functional  Q(\j;;x,t)  is  a  continuous  function 
over  AxRxR  where  the  admissible  set  A  is  given  the  weak  sequential 
topology. 

1  Proof:    The  result   follows  from   (3.10)    and      the      fact      Q(\|;;x,t)      is 

jointly    contluous   In  x  and   t. 

Theorem  3.6.  The  set  of  admissible  functions  A,  defined  by 
(2.17),  is  compact  in  the  weak  sequential  topology;  that  is,  every 
sequence   in  A   contains   a   subsequence   w-convergent   to  an  element    of   A. 

Proof:  By  LemnH  2,3,  the  total  variations  of  the  measures 
associated  with  A  are  uniformly  bounded.  It  follows  then  from  the 
Helly  selection  principle  that  every  sequence  of  associated  measures 
contains  a  subsequence  {\|^^(ri  )dri}  convergent  to  a  measure  da  in  the 
sense 


11m/    x(n)   ^]^(r])    dn    =   /    x('^)    do 
k-x" 


for   every   continuous  x*   We   conclude   from   (2.17)   that   for  all 
nonnegatlve  x 


0  <  /  X  (n  )  da  <  /  X  (n )  *  (n  )  dn 


This  implies  that  the  measure  do  satisfies 
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0  _<  da  _<  <|)  (n  )dri   , 

from  which  we  conclude  by  the  Radon-Nikodym  theorem  that  da  =  tj)  (n  )dn 
where  ^  satisfies  (2.17).  Thus  i|;.  is  w-convergent  to  \|)  e  A  and  the 
proof  of  Theorem  3.6  is  complete. 

Combining  the  compactness  of  A  asserted  by  Theorem  3.6  with  the 
continuity  of  Q  stated  in  Corollary  3.5,  it  will  follow  from  a  general 
topological  principle  that,  for  each  x  and  t,  Q(i);;x,t)  takes  on  a 
minimum  value  over  A.  This  gives  a  direct  argument  that  the  minimum 
problem  (2.16)  of  Theorem  2.2  has  a  solution. 

To  prove  uniqueness  of  the  solution  we  need  another  property  of 
the  operator  L. 

Theorem  3.7.  The  integral  operator  is  negative  definite  over  all 
i()  in  lMo.I]  for  which  Li()  is  in  L°°(R).   This  means 

(3.14)  -m>,^)  >  0 

for  all  \);  /  0  in  that  class. 

Proof:  It  is  convenient  to  extend  the  functions  \\) ,  originally 
defined  on  the  interval  [0,1],  to  all  R  by  first  setting  \})  (n  )  =  0  for 
n  >  1  and  second,  making  i|)  odd  by  setting  <|)(-ri)  =  -<|'(n).  This  allows 
us  to  rewrite  the  operator  L  defined  by  (2.22),  as  a  convolution 
operator. 


(3.15)  LiKn)    =-^     /    log(n-M)^  'i'(M)   dvi    =    -  £    *  ij;  (n  )    , 

ZTT 

— oo 

where 
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(3.16)  Jl(n)   =  "4-  log  'i^- 


Since   Li|;(ri)   is    odd,    we  have 


(3.17)  -(l.^^,,^)    =  1  <i*i|;,i);> 


where   the  brackets   on  the   right   denote   the  usual   duality   on  R, 

00 

<Xyi>'>  =    /  x(n)  (^di)  dn   . 

The  def initeness  of  convolution  operators  is  the  subject  of 
theorems  going  back  to  Caratheodory ,  Bochner  and  L.  Schwartz.  By 
(3.17),  Theorem  3.7  follows  from 

Lemma  3.8.  If  i|;  e  LHR)  is  odd  with  compact  support  and 
SL    *  i>  e    L"(R)  then 


(3.18)  <Z*^^,,^l>>    =  /  --L-  |il;(s)|2ds  , 

_  Is  I 


where  \(^(s)  is  the  Fourier  transform  given  by 


(3.19)        Us)   =  L^y  /  e-i^^  i|;(n)  dn  , 

)      -00 


(2Tr 


Proof:  Since  £  is  a  tempered  distribution,  for  all  ij;  in  the 
Schwartz  class  S  the  Fourier  transform  of  the  convolution  £*ij;  is 
/2tr   iU      By   the  Parseval   relation 

(3.20)  <li*^^),^l)>   =^2i^<ll,^>   =  /27  d,  li[i]^>   . 
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The  Fourier  transform  of  I    in  the  sense  of  the  theory  of  distributions 
is 


-^         1         1 

Jl(s)  = W  F.P. 


(27r) 


1/2      |s| 


where    the   F.P.      indicates   the   integral    is   a  Hadamard    finite   part, 
In   particular   if  x(0)    ~  0  then 


(3.21)  a.x>  =7^   /  ^x(s)  ds 


"  2 

We      apply      this     now     to  x(s)    =|i|^(s)|    ;    taking  \|;    to  be   an   odd   function, 

(3.19)   shows 


1(0)   =  L^     f  ,|;(n)  dn  =   0 

(2Tr)l/2   ^ 


and  allows  us  to  employ  (3.21)  in  conjunction  with  (3.20)  to  obtain 


(3.22)  <e*ip,i|;>  =  /  JL  i;i;(s)|2d5 

..     I  s  I 


for   all  odd  \|j    in  S  . 

Now    let  f    be   any   odd   compactly    supported  member  of  L   (R)   for  which 
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£*ij)   is   in  L"(R).   We  approximate  \)^  by  a  sequence  4;   of  functions  in  S 
defined  as  follows: 


(3.23)  ^n  =  Jn  *  * 


where   the  raollifier    j      is    given  by 


j^Cn  )    =  nj(nn) 


(3.24) 


j(n)   = ^rrexp  [-n^/2)    . 

(211)1/2 


Since  ^  is  odd  with  compact  support  and  the  j  are  even  elements  of  S, 
each  ^  is  an  odd  member  of  S  and  therefore  satisfies  (3.22).  By 
classical    properties    of    mollifiers   we  know  that    since  \|j    is    in  L    (R) 


(3.25)  lim  llij)j^-  i|)ll    1   =   0, 

rHK» 


and   since  £*>);    is    in  L°°(R) 


(3.26)  \iSi*i>J   „   =   II  j„*  £    *  il)ll    ^   <    U    *  i^ll   „ 


Using    (3.26)   and   the   fact  i  is   even,    see   (3.16),    we   get    the  estimate 

ia*(|)n,ij^j^>    -   a*ij;,ij;>|  =    |<£*i|;n+^*i|;,.|;n-^>| 

<  2IU*4;II       11,1.^    -   ,|;n     ,         , 
which  by   (3.25)    implies 
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(3.27)  lltn  <l*^^,yii^>   =  <S,*^^>,^>. 


Since   (3.22)    is    satisfied   for  ^^  we   deduce   from   (3.27)   that 


(3.28)  lira     /  JL   lil,    (s)|2ds   =    Ol*4,,i|<>   . 

n-HJO     -00 


However   frora    (3.24)  we   know 


^  1  9  9  \ 

Jn^^^    "  TTT  ^^   i~  s    /2n  J    , 

(2ii)^/2 


which  along  with  (3.23)  gives 

^^(s)    =  /27  jj^(s)  ii)(s) 

=  exp  (-  s2/2n^)  ^is)    , 

so   the   integral  on  the    left  hand   side   of   (3.28)    converges      to      integral 
in  (3.18),    completing   the  proof    of   both   the    lemma  and  Theorem  3.7. 
An    important    consequence   of   the   negative   definiteness   of   L   is: 
Theorem  3.9.        For     each      given     x     and    t,    the    functional  Q(\|);x,t) 
assumes    its  minimum   at    exactly    one   element,      denoted     ij;    (x,t),      of      the 
admissible    set   A. 

Proof:  Every  \|;  in  A  satisfies  0  _<  i);  _<  (})  and  (J)  is  in  L  [0,1]  so  we 
know  4)  is  also  in  L  [0,1].  By  Cbrollary  2.7,  Ijji  is  an  element  of  Cq(R) 
and  therefore  also  L°°(R).  We  can  then  apply  Theorem  3.7  to  the 
difference   of   functions    in   A   to  prove   that   Q(i|);x,t)    is    strictly      convex 
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over      the      set     A.    Since   A  is  a   convex   set,    Q(\|; )   assimes    its  miniraum  at 
exactly   one  admissible  4> . 

Theorem  3.10.       (a)   ^    (x,t)    depends    continuously   on     x,      t      in      the 
sense   of  w-convergence . 

(b)    Lif)    (x,t)    depends    continuously   on  x,t    in   the   maximum  norm. 

Proof :      Suppose    (x    ,t    )  ->■    (x,t).      Since,    according  to  Theorem  3.6, 

•k 

A  is  compact,  the  sequence  \|;  (x^^.t^)  has  a  cluster  point  ip*  in  A;  a 
subsequence  of  ij;  (^n'^-n^  w-converges  to  ij;*.  Since  by  Corollary  3.5, 
Q(i|j;x,t)    is   a      continuous      function     of     t,      and      since      by      definition 

■k  "k 

Q    (Xj^,t^)    =  Q(il)    (Xj^,tj^),Xj^,t^),    we    conclude   that    over   this   subsequence 


(3.29)  lim  Q    (x^^.t^^)    =  lim  Q(i|;    (x^^,  t^);  x^.t^,)    =  Q(t*;x,t). 

On  the  other  hand,  since  by  part  (b)  of  Theorem  3.1,  Q   is  a  continuous 
function  of  x  and  t 


lim  Q*(x^,tj^)    =  Q*(x,t). 

•k 

Comparing  this  with  (3.29)  we  see  Q  (x,t)  =  Q(4)*;x,t)  which  means  that 
t*  solves  the  minimum  problem  at  (x,t)  and  so,  by  the  uniqueness  of 
Theorem  3.9,  ij^^^  =  ij;  (x,t).  Since  i()  (x,t)  is  therefore  the  only  cluster 
point    of  <|;    (x    ,t    )   and  A  is   compact,   we  conclude 


•k  k 

w-lim  rp    (x^.tj^)    =  ^    (x,t)    ; 

this  proves  the  continuity  of  \li  .  The  continuity  of  lip  follows  since, 
by  Theorem  3.4,  L  is  sequentially  continuous  from  A  into  Cq(R)  and  the 
proof   is    complete  . 
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The  argument  given  in  Theorem  3.1  for  the  continuity  of  Q  (x,t) 
shows  that  Q  is  even  Lipschitz  continuous  in  x  and  t.  But  now  we  can 
show  more    is    true.  ' 

Theorem  3. 1 1.  Q  (x,t)  is  a  continuously  dif ferentiable  function 
of    X  and    t  with        ,  •  ...■■-'. 


3„Q*(x,t)    =       -  (n,i^*(x,t))      , 


(3.30) 


9^Q  (x,t)   =  - —  (n'^.tl^   (x,t)). 

I-  TT 


Proof:    We    deduce   from  formula    (3.1)    that 


(3.31)      Q(t;x,t)    -  Q(iJ;;y,s)    =  -  (n  ,i^)  (x-y)   -  i-  (n^,i|<)(t-s) 

IT  TT 


for   all  i|;  e    A.    Since  \|;      minimizes  Q,    we  have    the   two-sided   inequality 


Q(i|;*(x,t);x,t)-Q(i})*(x,t);y,s)  _<  Q*(x,t)    -Q*(y,s)      . 

<  Q(i|;*(y,s);x,t)    -  Q(i|;*(y  ,s)  ;y  ,s)      , 


which   when   combined  with   (3.31)    gives 


1  (n,^*(x,t))(x-y)   -I^  (n^.(;*(x,t))(t-s)  iQ*(x,t)   -  Q*(y,s) 

<  1  (n,i(^*(y,s))(x-y)  -—  (n^«|'*(y,s))(t-s). 

~  IT  H 
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Applylng   the   continuity   of  \()      in   the  w-topology  proved   in  Theorem     3.10 
to     both      sides      of      the  above   two-sided    inequality  yields   the   formulas 
(3.30)      and      shows      Q        is      continuously      dif f erentiable,      proving      the 
theorem. 

The   solution  of   a  minimum  problem   satisfies  variational   conditions. 
For  a  minimum  problem  with  constraints,    such  as   the  one  we  are  dealing 
with,    the  variational   conditions  take   the   form    (suppressing  x  and    t) 


(3.32)  -1  Q(^*  +ex)l  >   0 

de        "^  e  =0  - 


for   all   functions  x    such   that  \|j   +x   is    admissible.      Since 


Q(^4^X)    =  QC^')    +e   -  (a-L4<,x)    "  e^  -  (Lx.x). 


we    can  restate   (3.22)    so 


(3.32)'  (a  -  L^*,x)    >  0. 


Definition.        We     say      that      \j;      in     A      satisfies      the      variational 
condition  for  minimizing  Q(>|;)   in  A   if 


(3.33)  ^Q(^+ex)l      ^=(a-L^,x)2  0 


for   all  X    such   that 

(3.33)  '  i|;   +  X     e  A 
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Suppose  ip    satisfies   (3.33);   set 


X+  =     ■ 


-  \p     where      a-L\|;   >    0 


0        otherwise    > 


X-  = 


i(;     where     a-LiJj   <    0 


0  otherwise 


Clearly,   x+  and  x_  satisfy  (3.33)'.      By   construction,    (a    -  Li|;  ,x+ )  _<     0; 
this    is    the   opposite   of   (3.33),    which   implies    that 

.■.■.,~;       -      ;       ■    '  (a-L4),x+)    =0. 

This    implies    that 


(3.34) 


0     where      a-Lij;   >    0 


())     where      a-Li);   <    0 


It   is   easy   to  verify   that,  conversely,  if  ij;  satisfies  (3.34),  then 
(3.33)  holds  for  alli|;,  satisfying  (3.33)'. 


Theorem  3.12.   If  ij;  in  A  satisfies  the  variational  condition,  then 


"l^  =  t  • 
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Proof:    Since   x   =  "i**"^   satisfies   (3.33)',   \J;    satisfies   (3.33)  with  x 


*%, 


{a  -  L\p,   \l)      -  \\))    >  0   . 


Since  the  minimizing  i|)   also  satisfies  the  variational  condition,  and  x 


^~i>      satisfies  (3.33)',  it  satisfies  (3.32)': 


(a-Li|;*,i|;-ii)*)  >  0  . 


We  add  these  inequalities  and  obtain 


(L(iJ;*-iJ;),i)j*-((j)  >   0. 


Using   the   negative   definiteness   of  L  we    conclude   that  i|;    -rj;    =   0. 

The  importance  of  Theorem  3.12  lies  in  this:  in  order  to  show  that 
^  solves  the  minimum  problem  (2.16)  it  suffices  to  verify  that  it 
satisfies  the  variational  conditions  (3.34).  We  will  exploit  this  idea 
in  Sections   4    and   5. 

An  easy  direct  consequence  of  the  variational  conditions  (3.34) 
that   we  will  need   in   later   sections    is 

Theorem  3. 13.      For  all   x  and    t 

•5 

0     where     nx-4n    t  -  e,(n)   >  0 


(3.35) 


^    (x,t)    = 


where     rix-4n    t   -  9  _(ri  )    <  0. 
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Proof:   By   (2.42)   of   Corollary   2.7 


it 

a  <   a-L\);      <  a-L(|) 


Since      by     LemnH      2.6,   -Lij)  (n )    =e^(ri)    -e_(n)   and  by   definition   (1.22), 
a(n,x,t)    =  nx-4n    t   -  e^(ri),    we    see   the   inequality   can  be    rewritten  as 


nx-4n-^t  -  e^.(n)  _<  a-Li|j     <^  nx-An-^t  -  e_(n). 

Clearly   (3.35)   follows   from  this   inequality  and   the   variational 
conditions  (3.34);  and  the  theorem  is  proved. 

Titivated   by   the   variational  conditions,  we  partition  [0,l]xRxR 
into  three  disjoint  subsets: 

1°  =  {(n,x,t):  a  -  L;);*  =  0} 
(3.36)         I"^  =  {(n,x,t):  a  -  Li|;*  >  0} 

I"  =  {(n,x,t):  a  -  L.};*  <  0  }  . 


JL 

From  Theorem  3.10  we  see  that  a-L\|;  is  continuous  on  [0,l]xRxR  and  so 
the  sets  l"*"  and  I~  are  open  while  1  is  closed  and  separates  I  and  I  . 
We  denote  the  (x,t)  slices  of  I*^,  l"^,  and  I"  by  I°(x,t),  l"'"(x,t)  and 
I~(x,t)  respectively  where,  for  instance 

(3.37)         l'^(x,t)  ={ri:  (n.x.t)  e  1°}. 


Given      the   sets  I°(x,t),    I"'"(x,t)   and  I~(x,t),    one   can  determine  \|;    (x,t) 
by 

Theorem  3.14.      If   ij;  e    L^O,!]    and  Li};  e    L~(R)   such   that  \|;    satisfies 


and 
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CS.38)  LiJ;   =  a      on     I°(x,t)      , 


(3.39) 


\|)   =    0  on  I"'"(x,t) 


\j;    =  (J)   on  I    (x,  t) 


then  ij;    =  \()    (x,  t)  . 

if 

Proof:    That  i|;    (x,t)   satisfies   (3.38)   and    (3.39)    is   clear    from     the 

construction  (3.36)  of  the  sets  I  ,  I  ,  and  I~,  and  from  the 
variational  conditions  (3.34).  Now  if  41  satisfies  (3.38)  and  (3.39) 
then  the  difference  \\)q  =  ij;-\j;  (x,t)  satisfies  the  homogeneous  version  of 
(3.38)   and   (3.39): 


(3.40) 


L\|;q  =   0  on  I°(x,t)    , 


^        =   0  of f  l'^(x,t). 


This   implies    that    the   product  i|)q«L\|)q   =   0  for  all  t\  ,    so   that 

But,    according    to  Theorem  3.7,    L  is  negative   definite    so  we   conclude  iJ^q 
=  0,   proving    the    theorem. 

We     now     establish      the      last   major  property    of   the   operator  L.    We 
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extend  the  functions  ip  defined  on  [0,1]  to  all  of  R  by  setting  ij;  (n )  =  0 
for  n  >  1  and  taking  tj)  to  be  odd:  i|;(-ti)  =  -  i|^(n).  "Hie  set  of  all 
functions  obtained  by  extending  elements  of  LP[0,1]  (p  >^  1)  in  this 
fashion  we  denote  Lg(R).  Clearly  Lg(R)  is  just  the  set  of  all  odd 
lP(R)   functions   supported    in    [-1,1], 

Theorem  3.15.      If  i|;  e    Lg(R)   for   some   p,    1    <  p   <  »,    then  Li|)  e    Cq(R) 
and 

n 

(3.41)  L,|;(n)    =  /    H^(t)    dT    , 

0 

where  H   is    the  Hilbert    transform: 

00 

(3.42)  H,|;(n)    =Il^     /    -LiKm)    du. 

Proof:   Since   the  Hilbert   transform  takes  L^  functions  into  L^ 
functions  f or  1  <  p  <  »  such  that 


H^(t)  =LP-limi  /  ^"^   ^  i}/(w)  du. 


e^O  ""   -o  (Tni)^-^t2 


we  conclude 


/  H.|;(t)  dx  =  lim  /  i  /  I^IH -,|,(u)  dy  dT 

0  e->0  0  ^  ^  (t-u)^+c'^ 


The  order  of  integration  can  then  be  exchanged  by  Fubini's  theorem  and 
the  T  integration  can  be  carried  out  to  obtain 
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■n  00  9         7 

/    H^(T)dT    =   lim     /   i-log    ^"^^^   !^      ii,(y)    du 
0  e+0  -»    ^"^  P^^^ 


(3.43) 


=   lim       J-     /    logilinLl!^^(p)   dy 


Here   the   oddness    of  ^   was   used   in   the   last    step.      Since   for   fixed  r\ 


e.O^-  (n^)2^2        2.         ^^-^^ 


where   1/p  +    1/q   =   1,    we    see   from   (3.43)    that 


Tl  1 

(3.44)      /    H,j;(T)   dT    =^/    log  (11^)2  ^(y)    dp 

0  ^"^    0  ^^ 


which  by  the  definition  of  L,  (2.28),  proves  (3.41). 

The  continuity  of  Li|)  follows  directly  from  formula  (3.41).   It  can 

be  seen  from  the  definition  of  L  that  the  kernel  -—  log  f 1   tends  to 

2tt      n  -Hi 

zero  uniformly  over  the  support  of  ^    as  In  I  tends  to  infinity.   Thus  Lvp 
vanishes  at  infinity  and  the  theorem  is  established. 

Recall  now  the  relation  of  the  Hilbert  transform  to  analytic 
functions:  Let  RP  (1  <_p  <  »)  denote  the  space  of  all  functions  f(c) 
satisfying 

(1)  f(c)  is  analytic  in  the  upper  half  plane  Im(^ )  >  0. 

(2)  There  exists  a  constant  C  >  0  such  that 


(   /  If  (n+iT)|P  dn)^/P  <  C,  for  all  T  >  0. 
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Such  functions      have      boundary     values      in     LP(R)      In      the      sense      that 
LP-lim  f(n-H.T)    =  f(n)   exists.. 

The    following  result    is    classical    [7]: 

Theorem  3.16.  If  i|;  e  lP  for  some  p,  1  <  p  <  "» ,  and  ^  is  real, 
then  i|j  and  H^;  are  the  real  and  Imaginary  parts  of  the  boundary  values 
of   a   function   in  HP.      That    is 

^^)(r))      =     Re(f(n))       , 
(3.45) 

HiJ;(ti)    =  Ira[f(n)) 

for   some    f  e    HP. 

Using  this  result,  the  solution  of  equations  (3.38),  (3,39)  can  be 
reduced   to  a  Riemann-Hilbert   problem  in  potential   theory. 

Rather  than  determining  \j;  (x,t)  directly  we  shall  determine  the 
partial  derivatives,  8  ij;  (x,t)  and  3  ^^i];  (x,t).  We  will  later  make 
assumptions  that  imply  that  t\i  depends  dif f erentiably  on  x  and  t  in  the 
w-topology  over  the  ri-variable  with  derivatives  in  lP[0,1]  for  some  p, 
1  <  p  <  2.  This  means  that  for  every  x  i^  C[0,1]  the  function 
(x  j'l'  (x,t))  is  a  differentlable  function  of  x  and  t,  and  there  exist 
functions  '|'^(x,t)   and  )J;j^(x,t)    in  lP[0,1]    such    that 


*> 


(3.46)  ^  ^ 
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Whenever  the  support  of  x  lies  inside  the  open  sets  I  (x,t)  or 
I~(x,t)  we  see  by  (3.39)  that  the  left  hand  side  of  (3.46)  vanishes. 
Therefore 

(3.A7)  'I'*  =  'I't  =  0  off  I°(x,t). 

For  any  continuous  function  Xj  ^x  is  also  continuous  and  so,  since 
L  is  a  symmetric  operator,  we  deduce  that  (x  jL*!^  (x,t))  =  (Lx  ,t  (x,t)) 
is  a  dif f erentiable  function  of  x  and  t.  By  (3.46), 


^^ixM*)    =  (Lx.l'x)   =  (x.L^-x) 


(3.48) 


•k  "k  -k 


Let  Iq  be  the  interior  of  the  closed  set  I  given  by  (3.36)  in 
[0,l]xRxR  and  let  lQ(x,t)  denote  its  (x,t)  slice  defined  as  in  (3.37). 
Choose  the  support  of  x  in  (3.48)  to  lie  inside  the  open  set  lQ(x,t); 
then  by  (3.38),  14;   =  a  on  the  support  of  x^ 


(X.^*  )  =  (X.a)  • 


Using  the  definition  (1.23)  of  a  we  deduce  that 


(3.49)    3^(x,L4^*)  =  (x.n),  at(x,L%)  =  -^(x.n^). 
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Combinlng    (3.48)   and    (3.49)   we   find 


(3.50)  Lij;*  =  n    ,        14*   =  -  4n^     on     Io(x,t) 


By  Theorem  3.15  both  sides  of  (3.49)  are  continuous  functions  of  ri  ,  so 
the  equality  may  be  extended  to  lQ(x,t),  the  closure  of  lQ(x,t)  in 
[0,1]. 

Our  later  assumptions  will  imply  that  lQ(x,t)  and  1  (x,t)  differ 
by  at  most  a  finite  number  of  points.  Then  lQ(x,t)  and  I  (x,t)  will 
have  a  common  interior,  denoted  I(x,t),  on  which  (3.50)  holds.  Since 
I  (x,t)  and  I(x,t)  differ  by  a  set  of  measure  zero,  we  can  replace 
l'^(x,t)    by   I(x,t)    in   (3.47).      We  summarize: 

\l)      is    in  LP[0,1]    for   saie   p,    1   <p   <  °° ,    and   satisfies 


LT})^(x,t)    =  n   on     I(x,t), 


(3.51) 


ii;*(x,t)    =   0  off    I(x,t) 


Similarly  \l) ^   is   in  lP[0,1]    for   some   p,    1   <   p   <«,    and   satisfies 


f(  3 

L,j)      =  -4n        on     I(x,t)    , 


(3.52) 

i|)*   =     0        off      I(x,t) 
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Remark.  Arguing  as  in  Theorem  3.14,  it  follows  that  Tp  and  ^p ^  are 
the  unique  L^  solution  of  (3.51)  and  (3.52)  respectively. 

We  call  (3.51)  and  (3.52)  the  differentiated  variational 
conditions.  Note  that  the  inhomogeneous  terms  on  the  right  (i.e.  n 
and  -An  )  do  not  depend  on  the  initial  data  u(x). 

If  we  extend  I(x,t)  to  the  whole  real  line,  R,  by  reflection  about 
zero  then  (3.51)  and  (3.52)  hold  for  the  extended  functions.  We 
differentiate  the  top  equations  of  (3.51)  and  (3.52)  with  respect  to  n 
on  the  open  set  I(x,t);  using  (3.41)  we  obtain 


(3.53)  Hij;*  =   1     on  I(x,t)  , 


and 


(3.54)  Hil^*  =  -12n^  on  I(x,t)  . 


Denote   by  f    the    function   in  RP  whose   real  part    is  iJj    : 


(3.55)  Re(f)    =  4-! 


We  claim  that 


(3.56) 


Im(f)    =  1      on     I(x,t)    , 


Re(f )    =  0     off    I(x,t)    . 


The      first      relation      is    obtained  by  combining    (3.53)    and  Theorem  3.16; 
the   second   relation   is   from   (3.51). 

Similarly,    if   g  denotes    the  RP    function  whose   real   part    is  ^l>^: 
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(3.57)  Re(g)      =     ^* 


then 


(3.58) 


Im(g)   =  -  12  n^      on     I(x,t)    , 


Re(g)    =        0  off    I(x,t)    . 


We  shall   exploit   this    fact    in   the  next    two  sections      to      construct 

explicit        solutions        <|'^(x,t)        and        ^iiAxyt)  of      the      differentiated 

variational  conditions.     We  shall   then     verify  that      the      function     \|) 

satisfies   the  original  variational  condition. 

~      — 7~  — T  *  * 

We  next   show  how  to  express   u,    u     and  u     in  terras  of  i|)      and  ^^: 

Theorem  3.17. 


4  * 

u(x,t)    =  _  (n.^l'v^'^.t)) 


(3.59)  ^(x,t)    =^  (n,iP*(x,t))    =  -^  (r^^l|^*(x,t))  . 


^(x,t)   =  -  y^  (n^,i|'*(x,t))  . 


Proof:    We    use    formula    (3.30)    of   Theorem  3.11: 


9    0*(x,t)    =       _  (n,i|'*(x,t)) 


3,Q*(x,t)    =  -i^  (n\i^*(x,t))    . 
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We      use      formula    (3.46)    to  differentiate    these   with    respect    to  x  and    t. 
We    get 


9xxQ*('^''=)    =  -  (n,**(x,t)) 


(3.60)  8xtQ*(x,t)   =-  (n,ii'*(x,t))  =  -ii  (n^.|;*(x,t))   , 


9^^Q    (x,t)    =   - —  (r\^,^^(yi,t)). 


If  we  compare  (3.60)  with  definitions  (2.51),  (2.54)  and  (2.56)  for  u, 

— 2"  — 3" 

u     and  u     we    conclude  that    (3.59)   holds,    proving    the    theorem. 

We     conclude   this   section   by  proving    (2.29).      Using    the    definition 

(1.25)of  (})  ,    the   fact    that  L  and  —  commute,      and      relation      (3.41),      we 

dn 

ha  ve  - 


L<j)    =    -L_$    =    --l-L(t)=H$ 

dn  dn 


According      to  Theorem  3.16,    Lfi    is    the   imaginary   part    of    the   function   of 
class  hP  whose   real   part    is   $.      But    that    function   is 

00 

f(n)  =    /    [(-u(y)-n2)^^^  -  in]  dy  , 

-^10  •     .  . 

as   may   be   seen  from  (1.13). 

Using    (1.16)   and    (1.24)  we    see   that 

im(f )   =  e_(n)  -  9+(n)     . 
This   proves    (2.29). 
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4 .    The  Solution  until  Breaktline    . 

The  variational  problem  (2.16-2.18)  contains  x  and  t  as 
parameters.  We  first  investigate  the  case  t  =  0.  We  make  the 
assumption,  justified  a  posteriori,  that  the  set  I  consists  of  a  single 
interval 

(^•1)  I  =  (-6,6)  , 

where  B    is   a  dif ferentiable   function  of   x. 

We  claim  that  the  solution  of  the  Riemann-Hllbert  problem  posed  in 
(3.56)   is 


(A. 2)  f(0      =     1   - 


(62-^2)1/2     ' 


where  the  sign  of  the  square  root  is  chosen  so  that 

(4.3)  i(62-n2)l/2  >  q  for  n  >  6- 

Clearly  f(c )  is  analytic  in  the  upper  half  plane.   On  account  of  (4.3) 

(6  2-^2)1/2  _  -i^   for   Ic  I   large, 
which,  by  (4.2),  shows  that 


(4.4)  f(0   =  0(_i_) 

IC|2 


Since  f(?  )  has  singularities  of  the  form  (^-6)"^' 2  on  the  real  axis,  we 
see  from  (4.4)  that  for  1  £  p  <  2  there  exists  a  costant  C  >  0  such 
that 


(    /     |f(n+lT)|P)l/P   <  C 
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for  all  T  >  0.  This  shows  that  f  belongs  to  HP  for  1  <_  p  <  2. 
For  n  real,  it  follows  from  (4.2)  and  (4.3)  that 


(4.5) 


and 


Re(f)   = 


(e2-n2)l/2 


on  I 


off   I 


(4.6) 


Im(f)   = 


on  I 


1  - 


n2-B2)l/2 


off  I  , 


where  the  square  root  in  (4.6)  has  the  same  sign  as  r\  .  These  relations 
show  that  conditions  (3.56)  are  satisfied,  Re(f )  is  odd  with  support  in 
I,    and  Im(f)   is    even.     We  set    (4.5)    into    (3.55): 


(4.7) 


(b2^2)1/2 


for      In  I    < 


=     i 


for     In  I    > 


Now   using   fomula   (3.41)    of  Theorem  3.15,    with 


H,!^;;   =    Im(f) 
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given  by    (4.6),    we   have 


(A.  8)  L^; 


for     0   <  n    < 


[   n    -/ 


B    (y2_g2)l/2 


d\i      for     n   >  3 


Carrying   out    the   integral    in  (4.8)  we    obtain 


(4.9)  Ly\,l 


for     0   <  n    < 


,   n   -  (n^-B^)^/^     for    n   >  B      . 


The   value   of  0    is   so   far  undetermined;    for   that   we    turn   to   fornula 
(3.59), 


(4.10) 


-        4    .  *. 

"  =  -  (n,  'i'x^' 

IT  ^ 


We    expect    that 


(4.11) 


u(x,0)    =  u(x) 


Combining      (4.10)      with      (4.11)      and      expressing     ij;^  in   (4.10)   by    (4.7) 
yields 


u(x)   = 


(4.12) 


4    ,        *v 
-  (n,ip„) 

7t  '^ 


^/ 


^   0  (32-n2)l/2 


dn  =  -  B'    , 
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'1 

so  we    set  6  (x)    =/-u(x).      Since   -u(x)    >  n      between  x_(ti  )   and   x^Cn  )     and 
-u(x)    <  n^  outside,    see  Figure   1,    forraula   (4.7)    may   be   rewritten  as 


(4.13)  ^^U,0)      =    1 


^  for     x_(n)    <  X   <  x^(n  ) 

:-u(x)-n2)^/2 


for     X  _<  x_(n  )   or   X  >^  x^(ri ) , 


This    completes   the    determination   of   i);    (x,0). 

To  determine  ^  (x,0)  itself  we  integrate  'I'  (x,0)  as  given  by 
(4.13).  Theorem  3.13  states  that  \|;*(x,0)  =  0  where  nx  -  e^(n)  >  0;  in 
particular,  it  follows  that  for  fixed  n,  'I'  (x,0)  =  0  for  x  large 
enough.      Thus 


;j;*(x,0)    =  -/   ii;*(y,0)    dy    ,    i.e., 


(4.14) 


i>    {x,0)    =   ■ 


0 

x+(n ) 
/ 

X 


for     X  >^  x^(n  ) 

^  dy      for   x_(n)    <  x   <  x+(n  ) 

-u(y)-n2)^/2 


(-u(y 


f or   X  <  x_(n  ), 


In  (4.14)  we  have  used  the  definition  (1.25)  of  A  (n ) : 


(n)  = 


x^(n) 

/ 
x_(n) 


(-u(y)-n2)l/2 


dy 
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Clearly   0  <^  >|;    (x,0)   ^  <}]    so4)*(x,0)    Is    in   the  admissible   set  A. 

We  shall  now  show  that  4;  (x,0)  given  by  (4.14)  satisfies  the 
variational  conditions  (3.34),  thus  justifying  our  assumption  (4.1). 
Using    (4.13)   we    may    rewrite    formula    (4.9): 


(4.15) 


Li^x^'^''^^ 


for      x_(n)    <  X   <  X ,  (n  ) 


.  n-(n^  -Hi(x))^/2     for     X  <  x_(n)     or     x  >  x ,  (n  ) 


It    is    easy   to   see   from  Theorem  3.13   that 


so  by   Theorem  3.4  we    see 


w-lim     i|)    (x,0)      =     0   , 
x+<» 

lim     L\})*(x,0)      =      0   . 

X-Voo 


We    then   obtain  by    integrating    (4.15) 


-  ^,*U,0)      =     /    L4.*(y,0)   dy 


(4.16) 


,S 


/   n    -    (n^   +u(y))^/2  jy 


-  nx  +  e+(n) 


x_(n) 


-  nx  +  e+(n)  -  /         (n^  +  u(y))^/^  dy  for  X  <  x_(n) 


for  X  >  x^.(n  ) 

for   x_(ri  )  _<  X  ^  x+(n  ) 
1/2 


In    (4.16)   we    have   used   the   definition   (1.16)    of  9.(n): 


i+(n)   =nx+(n)  +      /       n   -  (n^-Ki(y))^^2  ^^^ 
x+(n  ) 


Since   a(n  ,x,0)    =nx  -  9,(ti),    we    see   from   (4.16)    that 


(4.17)   a-Li|j    =     • 
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/         (n^   +u(y))l/2   dy 
x+(n) 
0 

x_(ri  )    , 
-       /  (n2   +u(y))l/2   dy 


for  X  >  x^(n  ) 
for  X 


(n  )  _<  X  _<  x^(n  ) 


for  X  <  x_(n )   , 


which      clearly      shows      that     i|)    (x,0)      given     by      (4.14)      satisfies      the 
variational  conditions    (3.34). 

According  to  Theorem  3.12,  the  only  admissible  function  that 
satisfies  the    variational  conditions    is   ij;    .      Thus  we   have   proved 

Theorem  4. 1.  \\)  (x,0)  given  by  (4.14)  solves  the  variational 
problem   (2.16-2.18)    at   t   =  0. 

Formula  (4.14)shcws  that  i()  (x,0)  does  indeed  depend  dif ferentiably 
on  x  and  i(;^(x,0)  is  given  by  (4.13).  As  shown  in  Theorem  3.17  the 
differentiability  of  ij;  implies  that  (4.10)  holds;  comparing  this  with 
(4.13)   we    conclude   from  Theorem  2.12: 

Theorem  4.2. 

u(x,0)    =  weak  L^-lira  u(x,0;e)    =  u(x)    . 


e  +  0 


(4.18) 

This  result  was  expected  (see  (4.11))  and  provides  justification  for 
the  replacement  of  the  exact  scattering  data  by  the  asymptotic  data. 
In  Theorem  4 . 5  we    shall   sharpen  Theorem  4.2. 

We   now  begin  our  investigation   of   nonzero   values    of   t.    If      in      the 
KdV  equation   one  sets  e    =    0    , 
(^.19)  Uj.-   6uu^  =  0 

results.  In  the  introduction  we  have  pointed  out  that  there  is  a 
breaktime,  t^  ,  such  that  equation  (4.19)  has  a  smooth  solution  for  t  < 
t^  but  not  beyond.  Its  value  can  be  explicitly  determined  from  the 
initial  data   u(x) : 
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t.  =  [6  max  Uy(x) ]~  . 
°  X   ^ 


Denote  by  u(x,t)  the  solution  of  (4.19)  which  takes  on   the   prescribed 
initial  values  u(x).  .^ 

Theorem  4.3.   For  0  _<  t  <  t, 

(4.20)  u(x,t)  =  weak  L^-lira  u(x,t ;e )  =u(x,t) 

e->-0 

uniforTTily    in   t. 

Proof:  We  shall  reduce  this  to  Theorem  4.1.  As  u(x,t)  evolves 
according  to  (4.19)  it  will,  until  breaktime,  remain  of  the  same  class 
as  the  initial  data  u(x).  We  shall  consider  now  how  the  asmptotic 
scattering  data  corresponding  to  u(x,t)  change  with  t.  Define  x_j_(n  ,t) 
and   x_(n ,t)    as    the   two    roots   of 


(4.21)  u(x+,t)    +  n^   =   0 


with   x_(n,t)    <  x_^(n,t)    for  n    in  (0,1).      By   replacing  u(x)    by   u(x,t)   and 
x^.  (n )   by   x+(Ti,t)    in   formulas   (1.12)    and    (1.16)  we    define 

x+(ti  ,t) 

(4.22)  $(n,t)    =  /  [-u(y,t)-n^)^/2   dy     ^ 

x_(n,t) 
and 

oo 

(4.23)  e+(Ti,t)    =  nx+(n,t)    +  /         (n2-Ki(y,t))^/2   ^y. 

x+(n ,t) 


Lemma   4.4.      For  0_<t    <   t^^ 


(4.24)  x^.(n,t)    =  x+(n)   +  en^t   , 
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(4.25)  <&(n.t)  =*(n) 

(4.26)  e+(n,t)  =  n+(n)  +  4n^t. 

Proof:  Differentiating  (4.21)  with  respect  to  t  and  using   (4.19) 
we  obtain 


9  t-x+  =  -  —  =  -6u  =  6n  ■ 


Formula   (4.24)   new   follows   upon  integrating  this  and  using  the  fact 

'St('^>'^)  ~  x+(ti).   Differentiating  (4.22)  with  respect  to   t   and  using 
(4.19)  gives 


9tHn.t)  =-|  /       '  dy 

X   (.  u-ri  ; 


x+ 


=  -3  / 


X 


7:^^'^'' 


(2u  -  4n2)(-  u  -  n^)^/^l   =  0 


by      (4.21).        This      proves      (4.25).        Preforming      the    same   procedure  on 
(4.23)    shows 


te+(n.t)   =  -y    /  ^\       dy 


"+ 


"    ~    ^       ^  .     2    "^l/2    "X    ^^ 

x^     (n^+u)i/^ 
=   (4n2   -   2u)(n2  +  u)l/2 1"      =   4^3 
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which  after    integrating  yields    (4.26)   and   establishes   Lemma   4.4, 
From   (4.25)    of   Lemma    4.4  we    see   that 


<ti(n)    =  -  9^$(n,t)      , 


and  using  (4.26)  in  the  definition  (1.23)  of  a,  that 

a(n,x,t)   =  nx  -  e^.(Ti  ,t) . 

Thus    the   variational   problem    (2.16-2.18)   for  0  _<  t    <   t.     is    the   same      as 
the  variational   problem   for   t   =   0  except   that    u(x)    has   been   replaced  by 
u(x,t).      By  Theorem  4.1,    its  solution    is    given  by   formula      (4.14)      with 
the   same    replacements    and  Theorem  4.3   follows   from  Theorem  4.2. 
The  weak    limit   in   formula    (4.20)    is   in    fact   a    strong  limit. 


Theorem  4.5.      For  0  <_  tg  <    t^, 


(4.27)  u(x,t)    =  L^-lim  u(x,t ;e  )    =u(x,t) 


uniformly    in   t. 

Remark.      In    particular   this   proves  Theorem   1.2. 
Proof:   LeranH    2.13  gives 
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00  oo 

lim  /  u^(x,t,e)  =  /  u^(x)  dx   , 

G-vO  -»  -« 

uniformly   in  t.   Since  /   u  dx  is  a  conserved  quantity  for  equation 
(4.19), 


(4.28)     lim  f  u2(x,t;e)  dx  =  /  u2(x,t)  dx  , 
e*0  "^  ^ 


uniformly  in  t  for  0  <^  t  <  tu . 
By  Theorem  4.3 


(4.29)     lim  /  w(x)u(x,t;e)  dx  =  /  w(x)u(x,t)  dx  , 
e+  0  -»  -«> 


uniformly  in  t  for  0  <^  t  <  t^^  and  over  w  in  any  compact  subset  of 
L^(R).  But  it  is  easily  seen  that  for  0  ^  t  <^  t,  ,  u(x,t)  is  a 
continuous  function  of  t  into  L'^(R)  of  the  x  variable.  Thus  the  set 
{u(x,t):  0  <.  t  <_  t^}  is  a  compact  subset  of  L  (R);  setting  w(x)  = 
u(x,t)   in    (4.29)    we   conclude   that 


(4.30)  lim     /    u(x,t)u(x,t;e )   dx     =     /   u^(x,t)    dx      , 

£+0      -OO  00 


uniformly    in   t    for  0  _<  t    <   t.  . 
Consider  the   identity 

00  oo 

/    (u(x,t;e  )    -  u(x,t)]      dx  =     /    u    (x,t  ;e  )   dx 

— oo  — oo 

00 

-    1     I    u(xjt)   u(x,t;G)   dx 

-00 

oo 

+  /    u^(x,t)    dx    . 
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By   (4.28)   and    (4.30)    the  righthand   side   tends   to     zero     uniformly      over 
0  _<   t    <   t^^  as   e    tends    to  zero,    proving   the   theorem. 

In  the  proof  of  Theorem  4.3  we  relied  on  the  knowledge  that  u(x,t) 
=  u(x,t)  satisfies  (4.19)  for  0  £  t  <  t^.  We  show  now  how  (4.19)  can 
be    deduced   directly   from   the    variational  conditions. 

Assume  that  the  set  I  consists  of  a  single  interval,  see  (4.1), 
where  g  is  a  dif f erentiable  function  of  x  and  t.  This  assumption  will 
be  justified  a  posteriori.  We  claim  that  the  solution  of  the 
Rieraann-Hilbert   problem  posed    in   (3.58)   of   Theorem  3.18   is 


(4.31) 


g(c)=-12,2i,i2^^^     , 


where  the  sign  of  the  square  root  is  again  chosen  to  satisfy  (4.3). 
The  same  arguments  that  applied  to  f(^)  given  by  (4.2)  can  now  be 
applied   to  g(c )   to   show  g   belongs    to  HP   for    1   <   p    <  2,   and    for  n    real 


(4.32) 


and 


Re(g)     = 


12n^-6g^ 
(B2-n2)l/2 


on     I 


off      I    , 


(4.33) 


Im(g)      =  - 


-   12  n' 


-  12  n    + 


3        2 

12T1  -ee'^n 


on     I 


off      I 


-74- 


As  before,  the  square  root  in  (4.33)  has  the  same  sign  as  n •  These 
relations  show  that  conditions  (3.58)  are  satisfied,  Re(g)  is  odd,  with 
support  in  I,  and  Im(g)  is  even. 

We  set  (4.32)  into  (3.57)  to  obtain 


(4.34) 


12n^-6e^n 
(b2^2)1/2 


for   In  I  < 


for  In  I  > 


Now   using   fonmla   (3.41)    of  Theorem  3.15  with 


Hi|^^   =   Im(g) 


given  by    (4.33)   we   have 


(4.35) 


4  n 


for     0  <  n    < 


^  '^      "•"  J    — ^ — ^    ,  ,^   dy    for  n    > 


3    (p2-32)l/2 


Carrying  out  the  integral  in  (4.35)  we  obtain 
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-  4  n- 


(4.36)  h^^>] 


for     0  <  n    <  B 


-  4  n^+  (4n2+2e2)(n2_g2)l/2       f^^     ^   >  g 


"k  "k 

Now  \\)      given  by   (4.7)   and  \|jj.    given     by      (4.34)      must      satisfy      the 
conpatibility   condition  9  ^\(;      =  Sv'^f*      "^i^    implies 


0  =^(-VtT7t)   ^3x('^'"fwJ 
(^2^2)1/2  ^   (g2^2)l/2 


(4.37) 


nB 


(B2-n2)3/2 


(^t  +  ^e^ej  , 


from  which  we    conclude 


(4.38) 


Pt  +  6  B'^Bx  =  °  • 


Expressing  i|;      in   fomula    (4.10)   by    (4.7),    we    obtain    (see  4.12) 


(4.39) 


4  * 

u(x,t)   =  -  (n ,^^)   = 


=  -  B^Cx.t) 


Multiplying    (4.38)    by  -2B    and  using    (4.39)  we    deduce 


(4.40) 


"t  -   ^"\  =  0 


this  is  relation  (4.19),  valid  for  t  <   t^. 

The  above   procedure  will  be  generalized  in  the  next  section  to 
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obtain   the   time   dependence   of  solutions    of   the   variational  problem     for 


t   >    t^. 


We  close  this  section  with  an  observation.   Set  \\)      given  by  (A, 7) 


and  ij;*  given  by  (4.34)  into  (3.59)  of  Theorem  3.17: 


uV.O=--<n3,*;,=^;_i_.,, 


u^(x,t)   =  -^  (n^ift)  =  "iz-  /        '     ri/9^^ 


A  brief   calculation  shows    that 

u^x.t)    =       B'^(x,t)    =u^(x,t)    , 
(4.41) 

u^x.t)    =   -6^(x,t)    =u^(x,t)    . 


—    — y         — T 

From    the    definitions    (2.51),    (2.54)    and    (2.56)    for  u,    u     and  u     we  must 
have 


(4.42) 


u^-    3  ^j^  =  0 


"  t   -   ^  "   X  =  0 


Using  (4.41)  to  express  u  and  u  in  terms  of  u,  we  see  each  of  (4.42) 
will  hold  if  and  only  if  (4.40)  holds,  showing  again  the  validity  of 
(4.40). 
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5.   The   Solution   beyond   Breaktlme. 

In  this  section  we  describe  the  solution  of  the  variational 
problem  (2.16-2.18)  when  t  is  larger  than  the  breaktlme,  t,  .  We  make 
the  assumption  that  the  set  1  consists  of  a  finite  union  of  disjoint 
intervals,    say 

(5.1)  I  =   (62n.e2n-l)        •••        (^2.31)     , 
where 

(5.2)  e2n  ^  ^2n-l  <  ...  <  ^2  ^  ^1  • 
Since  I  is  symmetric  about  n  =  0,  we  see  that 

^Zn-k+l  "  ~  \  ' 

and  that  the  set  1  Is  completely  determined  by  its  set  of  positive  end 
points 

(5.3)  B  =   {Bk>k=l    • 

We  will  work  as  if  n  is  odd  since,  as  we  will  later  see,  this  will  be 
the  case  except  when  the  initial  function  u(x)  takes  on  the  value  zero; 
in  that  case  we  make  n  odd  by  setting  6   =  0. 

As  first  step  we  shall  construct  all  possible   solutions   of   the 

"k  "k 

Riemann-Hilbert  problems  (3.56)  and  (3.58)  that  i|/  and  \l) ^  satisfy.  The 
solution  will  be  expressed    in   terms   of   a   function  K(c,)    defined   by 

n 

(5.4)  r2(0    =    I     I    (Bk^   -  C^)    ; 

k=l 

we    choose   its    sign   so   that 
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(5.5)         iR(n)  >  0  for  n  > 


Lemma  5.1.   All  functions  of  the  form 


P(c) 
(5.6)  f(c)  =   i  -  ^^^ 


R(0 


(5.7)  g(0  =  -  12  c^i  +  ^^^^ 


R(C)  ' 


are  solutions  of  the  Riemann-Hilbert  problems  (3.56)  and  (3.58) 
respectively,  where  P(5 )  and  Q(c)  are  odd  polynomials  with  real 
coefficients 

(5.8)  P(0    =  C""  +  a2?""^  +  ...  +  a^_2C 

(5.9)  Q(0    =   12^"+2   _   g^j-   e^2)^n  ^^^^n-2  +_  ^+  ^^_^^    ^ 

Remark.   It  is  not  hard  to  show  that  these  are  all  solutions  [9]. 
Proof:  Clearly  f(c )  and  g(c )  are  analytic  in  the  upper  half  plane. 
By  (5.4),  (5.5) 


(5.10)         iR(o  ~    ^""-jil    ej)?""2 +0(IU""^) 

k=l 


for  1^1  large,  which,  along  with  (5.6-5.9),  shows  that 


(5.11)       f(o=o(_!_),     g(o=o(-i_] 

Id  I?  r 
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Slnce  f(c)  and  g(c  )  have  singularities  of  the  form  (?-6)^'^  on  the  real 
axis,  we  see  from  (5.11)  that  f  and  g  belong  to  HP  for  1  <^  p  <  2, 

For  n  real,  the  function  R(n)  defined  by  (5.4)  is  real  and  even  on 
I,  imaginary  and  odd  off  I,  where  I  is  defined  by  (5.1).  Since  P(n  ) 
and  Q(n)  are  real  and  odd  for  n  real,  it  follows  from  (5.6)  and  (5.7) 
that: 


(5.12) 


Re(f)  = 


P(n) 
R(n) 


on  I 


off  I  , 


and 


(5.13) 


lm(f)  = 


on  I 


1+1  —7 — r-  or r  L    , 

RTnT 


while 


Q(n) 


on  I 


(5.14) 


Re(g)  = 


0   off  I  , 


and 
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(5.15) 


Ini(g)    = 


12   n 


on      I 


-n.^-.|^   -" 


These  relations   show  that    f(c)  and   g(c )   satisfy   conditions      (3.56)      and 
(3.58)   as   asserted   in  Lemma   5.1. 

Using   (3.55)   and    (3.57)  we    set 


(5.16) 


fx  = 


(5.17) 


-tg^l 


We  now  show  how  to  choose  the  coefficients  a^  and  y  •  so  that  the 
differentiated  variational  conditions,  (3.51)  and  (3.52),  are 
satisfied.   Using  formula  (3.41)  with 


given  by  (5.13)  and  with 


given  by  (5.15),  we  obtain 


K 


Im(f  ) 


Hij;^  =  Im(g) 


(5.18) 


M>: 


P(u). 


- 1  -(w^) '" 


and 


(5.19) 


L^,  =  -4,3  +/  Im[^)  dp 
0    ^^ ' 


-81- 
If    the   differentiated  variational  conditions,    (3.51)   and    (3.52)   are     to 
be   satisfied,   the   integrals   on   the   righthand   sides    of    (5.18)   and    (5.19) 
mist   vanish  on  I.    Since    the   integrands  are   even   functions   which     vanish 
on  I ,    we    see  it    is    both   necessary  and    sufficient   to  require 

62k 

(5.20)  /       ^d„  =   0     , 

P2k+1 
e2k 


Q(n) 
2k+l 


(5.21)  /       H^dn-0 


for  k  =  1,2,...,——  .  Recalling  the  form  of  P  and  Q  given  in  (5.8)  and 
(5.9),  we  see  that  (5.20)  is  a  system  of  (n-l)/2  inhomogeneous  linear 
equations  for  the  (n-l)/2  unknowns  012,04,  ...  ,aj^_l  and  (5.21)  is  a 
system   of      (n-l)/2      inhomogeneous     linear      equations      for     the      (n-l)/2 

unknowns  Y  2«Y4' '  ' ''^n-l    * 

Lemma  5.2.  The  systems  of  equations  (5.20)  and  (5.21)  have  unique 
solutions. 

First  Proof:  If  not,  then  according  to  the  alternative  principle 
of  linear  algebra,  the  common  corresponding  homogeneous  system  would 
have  a  nontrivial  solution.  But  then  there  would  be  a  nontrivial  odd 
polynomial  N(n )  with  real  coefficients  and  of  degree  at  most  n-2  that 
satisfies 


e2k 


N(n) 


(5.22)  /       ^dn  =  0 

ft  ^^^^ 

P2k+1 


for  k  =    1,2,...,  _— -  .      This   implies 
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"^0   =  Re  Ir^J 
satisfies 

Li|)Q  =  0      on  I 
\|;q  =  0      off  I  . 

According  to  the  uniqueness  argument  given  for  (3.40),  such  a  iJ;q  is 
zero,  and  so  N(ri  )  is  zero,  a  contradiction. 

Second  Proof:  If  not,  then  let  N(n)  satisfy  (5.22)  as  above.  But 
this  contradicts  the  following  lemma. 

Lemma  5.3.  Let  N(n )  be  any  nontrivial  odd  polynomial  with  real 
coefficients  that  satisfies  (5.22).  Then  N(n)  nust  have  at  least  one 
root  in  each  of  the  (n-l)/2  intervals  (B2k+i,e2k)  and  the  degree  of 
N(n)  Tiust  be  at  least  n. 

Proof:  Since  R(ri)  does  not  change  sign  in  (B2k+l»^2k^»  thenN(n) 
mist  if  (5.22)  is  to  hold.  Thus  N(ri)  has  at  least  one  positive  root  in 
each  of  the  (n-l)/2  intervals  (3 2k+l  »P2k^  •  Since  N  is  odd  it  has  an 
equal  number  of  negative  roots  and  a  root  at  zero,  altogether  a  minimum 
of  n  roots,  proving  the  lemma. 

According  to  (5.20),  the  polynomial  P(n)  satisfies  hypothesis 
(5.22)  of  Lemma  5.3.   Since  P  is  of  degree  n,  we  conclude 

Corollary  5.4.   P(ri)  has  a   zero   inside   each   of   the   intervals 

^^2k+l'^2k^'  ^  ~  l,...,n-l  and  at  ri  =  0,  and  nowhere  else. 

According   to  Lemma   5.2,  we  can  solve  (5.20)  and  (5.21)  uniquely 

for  a-  and  Y -i  ,  j  =  2,4,...,n-l  as  functions  of  the  points  &^  in  B,  see 
(5.3).  Observe  that  the  coefficients  and  the  inhomogeneous  terras 
appearing  in   the   linear   systems   (5.20)   and   (5.21)   are   complete 
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hyperelliptic  integrals.   Thus  the  a.  and  y ^  are  ratios  of  determinants 
whose  entries  are  such  integrals. 

We  have  now  proved 

Th e ore m  5.5.  Given  the  set  B,  \j)  as  determined  by  (5.16)  and 
(5.20)  is  the  solution  of  the  differential  variational  conditions 
(3.51)  and  4)  as  determined  by  (5.17)  and  (5.21)  is  the  solution  of  the 
differentiated  variational  condition  (3.52). 

The  functions  \|;  and  ijj^  satisfy  the  compatibility  condition  9^.1};^  = 
9  ij^j.;  or  equivalently 


rs  23->        a  P(^)  +  a  Q(^)  _  n 


We    shall  see  what    conditions   this    imposes    on   the   set   B, 
Rewrite   the   left   side    (5.23)   as 


(5.24)    3,   |<|>.  3,  «i> -3 


We  claim  that  the  numerator  on  the  right  is  an  odd  polynomial  of  degree 
at  most  3n-2.  For  from  the  definitions  (5.4),  (5.8)  and  (5.9)  we  see 
that  for  large  lc  I 

r2(c)  =  (-C^)"  +  [  I     3^)  (-C2)""^+0(IU2''"^)  , 

k=l 


?U)     =  ^"  +  o(ici""2) 


Q(C)   =   12  C"''^  -  6(  I     e2)  5n  +  o(lcl"-2) 

k=l 


so  that 
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n 

(5.25)  r2(3^p  +  3^Q)    =    (-^2)n  9^(_6     ^      g2)^n  +  o(|^  |3n-2)     ^ 

k=l 

and 

(5.26)  1  (P3tR2+  09^2)    .   1    (i2^n+2)   3^^      ^      gj)  (-^  2)""!   +  0(  I?  I  3n-2) , 
^  k=l 


Subtracting    (5.26)   from   (5.25)   we    find   the   leading   terms    on     the      right 
cancel,    demonstrating   the   claim. 

Since  R^  =  0  at  ^  =  B]^  ,  the  numerator  on  the  right  of  (5.24) 
vanishes  at  ^  =  g,  if  and  only  if  its  second  term  vanishes  at  ^  =  &\^  • 
That   second   terra  may   be  written 


n 
i  (P9^r2  +  Q9^Jl2)    =     ^      [Bj^[p(09t6k  +  QU)9xek)      "      (Bj-C^)]    • 

k=l  ji^k 


Clearly  this  vanishes  at  ^  =  g.  if  and  only  if  Q^   satisfies 
(5.27)     P(Bj^)  3j.B^  +<)(Bk)9xPk  =  ^  • 


This   time   evolution  of   the  B^   allows   us  to  factor  R  (c )  from  the 
numerator  of  (5.24)  and  write 


(5.28)  3,1^+3,^  =  1^   , 

t  R(C)    ^  R(C)   R(C)   ' 


where  N(^)   is    the   odd  polynomial   of   degree  at  most   n-2   given  by 


"  P(c)-P(Bk)  Q(c)-Q(6k) 

N(c)    =3,P(0    +9,Q(0    +     I  X ^-Bk3tek+ o 2~^k9xek 

k=l  C27B^2  ^2^    2 
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In    fact,   we  will  show  N(^  )    =  0. 

For  e    sufficiently    small  we      employ      the     fundamental     theorem     of 
calculus   to  obtain 

32k^  B2k-^  ,    , 

P  2k+l^  P  2k+l^ 

^2k-^ 


J         t  R(n)         X  R(n) 


(5.29)  ^2k+l'^^ 


^P(g2k^)  3te2k+  Q(e2k-E)  M2k 


» 


R(e2k-e> 

_  P(g2k+l-fe)  3tg2k+l+  Q(g2k+l->^)  9xg2k+l 
R(e2k+1-^) 

for  k  =  1 ,2, . . . ,(n-l) /2.  As  e  goes  to  zero,  the  denominators  of  the 
second  and  third  terms  on  the  right  of  (5.29)  will,  by  (5.4),  vanish 
like  e  '  ,  but,  by  (5.27),  the  numerators  of  the  same  terms  are  0(e); 
thus  the  limit  of  those  terms  as  e  -»■  0  is  zero.  Using  (5.28)  in  the 
first    term,    we    then  obtain 

B2k^  ^2\-^ 

(5.30)      lim[3,       /  l^n-^8,       /  §^  dn] 

e^O         e2k+l-te  ^2V.+\^ 

^2k 

P2k+1 

uniformly  on  compact  subsets  of  x  and  t.  On  the  other  and,  by  (5.20) 
and    (5.21), 
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lira         /         I^  dn   =   0   , 


(5.31) 


lira  /  §^  dn    =    0    . 


uniformly    on  compact   subsets    of   x  and   t.      Since      the      divergence      is      a 
closed   operator,   we  conclude   from    (5.30)    and    (5.31)    that 


^2k 
r        N(n) 

R  R(n) 


dn  =   0  . 


Using  Leraraa   5.3,    we    then   conclude  N(n)    =   0.    Thus  we   have   shown 

Theorem  5.6.         If    the  3j^   in  B   satisfy    (5.27),    then  ij;^  and  ip^   given 

by   (5.16)   and    (5.17)   are   compatible. 

We  deduce   from  Corollary   5.4  that 

P(Bj^)    ^  0    ,        k  =    l,...,n    . 
Dividing    (5.27)    by   P(3j^.)    gives 

(5.32)  9^.6  +  V(e;B)    3^B     for    B      in  B, 

where 

Q(n) 


v(n;B)     = 


P(n) 


indicates    the    dependence   of   V   on  the    set    B   through   the   dependence   of  a  • 
and  Y  .   on  B . 

Formula      (5.32)      is     a      coupled      system     of   n   partial  differential 
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equations  where  each  6  in  B  is  a  Rieraann  invariant  with  corresponding 
characteristic  velocity  V(B;B).  The  set  B  evolves  according  to  (5.32) 
so  long  as  the  solutions  remain  regular  and  n  does  not  change.  We  have 
not  yet  said  how  the  set  B,  and  its  cardinality,  change  when  a 
sigularity  is  encountered,  nor  have  we  addressed  how  B(x,t)  is  related 
to  the  initial  data  u(x).  We  shall  deal  with  these  issues  by 
introducing  a  pair  of  functions  x, (n ,t)  and  x_(n,t),  0  ^n  ^  1,  already 
defined  in  Section  4  for  t  <  t,  but  here  extended  to  all  t  ;  we  shall 
construct    the   set    B   in  terms   of    these    functions. 

We  start    by   defining   the   notion   of   crossing: 

Definition.  A  pair  of  functions  x,(«),  x_(« )  cross  the  value  x  at 
n  if  the  union  of  the  images  of  any  neighborhood  of  ri  under  x^.(»  )  and 
x_(* )  is  a  neighborhood  of  x.  Such  an  n  is  called  a  crossing  point  of 
the   pair   at    x. 

Note  that  if  0  <  n  <  1  and  x,.  (n )  =  x,  then  n  is  a  crossing  point 
at   X  unless   x^ (•  )  has   a    local   extreraura   at  n. 

For  a  given  pair  of  functions  x^(n),  x_(n )  and  a  given  value  x  we 
define    the   set   B(x)    by 

B(x)    =   {ri  :   n    is   a   crossing  point   of   x^(«),   x_(«  )  at    x}    . 

Note      that      if      u(x)      is      our   initial   function  and   Xf  (n )   are   defined  by 
(1.11),    then  B(x)    consists   of    the   single   point  3 (x)    defined  by    (4.13): 


6(x)    =  [-u(x))l/2_ 
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Suppose  the  functions  x^.  =  x+(«,t)  depend  on   the  parameter   t  ; 
then  so  does  the  set  B  defined  above: 

(5.33)  B(x,t)  =  In;  n  is  a  crossing  point  of  x^(.,t),  x_(«  ,t)  at  x} 

We  formilate  the  following  initial  value  problem  for  the  functions 
x+(ri  ,t)  ,  x_(n  ,t)  : 

dx^ 

— —  =  v(n;B(x^,t))  ,  x^(n,0)  =x+(n) 

(5.34)  i!^=  v(n;B(x_,t))  ,  x_(ti  ,0)  =x_(n), 

at 

where  x^(n  )  and  x_(ri  )  are  defined  by  (1.11).  Note  that  the  set  B(x,t) 
depends  on  the  functions  x,  ,x_  .  As  long  as  that  dependence  is 
Lipschitz  continuous  in  a  suitable  norm  for  x_j.  ,  x_  ,  (5.34)  has  a 
unique  solution.  We  assume  that  (5.34)  has  unique  solutions  such  that 

(a)  As   functions   of  n   and   t,   x^  and  x_  re  C   in  (0,1)  x  R  and 
continuous  in  (0,1]  x  R. 

(b)  The  limits  lim  x+(n,t)  exist,  possibly  as  ±  " .   If  finite  we  call 

n+0 
it   the   boundary   value  of  x^.  at  n  =  0,  denote  it  by  x+(0,t)  and 

assume  it  is  C*^  in  t. 

(c)  The  number  of  critical  points  of  x, (• ,t)  and   x_(«  ,t)   is   finite 
for  all  t. 

(d)  If   for   some  t\    ,  x^(ri )  =  x_(n),  then  x^(n,t)  =  x_(ri  ,t)     for  all 

t. 

Note  that  (c)  implies  that  the  cardinality  of  B(x,t)  is  finite. 
Since  x  (1)  =  x_(l) ,  it  follows  that  x  (l,t)  =  x_(l,t)  ;  we  denote  this 
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comiTDn  value  by  XQ(t).   From  (d)  and  the  fact  that  x_(ri  )  <  x+(n  )  for  0 
<  n  <  1  we  conclude  that  x_(n,t)  <  x^(ri  ,t)  for  all  t  as  well. 

Note  that  according  to  the  definition  of  crossing,  n  =  1  is  a 
crossing  point  at  x  =  XQ(t)  iff  for  n  near  1 

x_(ri,t)  <  XgCt)  <  x+(n,t). 

Lemma  5.7.  Suppose  that  x  is  not  a  critical  value  of  x^(«  ,t)  or 
x_(.  ,t),  nor  a  boundary  value  of  x+(«,t)  (i.e.,  x  ^  3C^(0,t)  or  x  5^ 
XQ(t))   then   for  all  x,t    in  a   neighborhood   of    x,t 

(a)  The   cardinality   of    B(x,t)    is   constant 

(b)  B(x,t)    is  a   dif ferentiable   function   of   x,t. 

Proof:  It  follows  from  our  assumptions  that  we  can  find  a 
neighborhood  of  x,t  such  that  9^x+  are  bounded  away  from  zero  and  x  is 
bounded  away  from  the  boundary  values  of  x^.(»,t).  By  the  implicit 
function  theorem,  there  exists  a  (maybe  smaller)  neighborhood  N  of  x,t 
such  that  every  solution  of  x  =  x^(n,t)  is  given  by  n  =  6(x,t)  where 
B(*,*)  is  a  c'^  function  in  N.  The  conclusions  of  Lemma  5.7  readily 
follows. 

Now   construct  I(x,t)    fromB(x,t)    according   to   (5.1). 

Lemma   5.8.      For   any   x, t   and   any  n    in  (0,1) 

(5.35)  x_(n,t)    <  X  <  x+(n,t)    =>  n  e   I(x,t)     , 


(5.36)       X  <  x_(n ,t)     or     x  >  x+(n ,t)    =>  n  ^  I(x,t)   . 

Proof:  By  (5.1)  if  n  is  not  in  B(x,t)  then  n  is  in  I(x,t)   if   and 
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only      if    the  number  of   elements  of   B(x,t)    in    [ri,l]    is  odd.      If   x  ?*  ^o^^^ 
then  by    (5.33)    the  number  of    elements  of   B(x,t)    in    [ri,l]    is   equal    to   the 
number  of   times   that  x    (p,t)    or  x_(M,t)    crosses   the  value  x   for   y    in    [n,l] 
Consider   the   case  x_(ti  ,t)    <  x   <  x^.(n,t);    if   x  <  XQ(t)      =     x+(l,t)      then 
x^(li,t)      mast      cross      the   value   x  at   an  even  number   of  y    in   [n,l]   while 
x_(p,t)    must   cross    the   value   x  at   an  odd  number  of  y    in      [n,!].        Thus, 
the      total     number  of    crossings    beings    odd  we    conclude  n   e    I(x,t).      The 
argument    is   essentially    the   same    for  all   the   cases,    provided   x  j^   XpCt). 
In    that   case   one  imst   consider   the   behavior  of    x^.  and  x_  near  p    =    1  but 
it    is    .iust   as   easy. 

Now   we  use    formulas    (5.16)   and    (5.17)      to     construct      'i^jj(x,t)      and 
i|;j^(x,t)   from  B(x,t). 

Lemma  5.9.  If  x  is  not  a  critical  or  boundary  value  of  x^(«  ,t)  or 
x_(«  ,t)    then  in  a   neighborhood   of    x,t 

(a)  ijj    (x,t)   and  \|j^.(x,t)   are   dif f erentiable   functions    of    x,t. 

(b)  The   compatibility   condition 

^t*x  =  ^x^t 
is  satisfied. 

Proof:  The  coefficients  a.  and  y  ■?  of  P  and  Q  are  functions  of 
3i,...,3j^  ,  expressible  as  rational  functions  of  complete  hyperelliptic 
integrals.  Therefore  they  are  analytic  functions  of  g  j^ , . . .  ,g  except 
if  two  g,  are  equal.  It  follows  from  Lemma  5.7  that  this  does  not 
happen  under  the  hypothesis  of  Lemma  5.9.  It  does  not  follows  from 
Lemma  5.7  that  the  3^  are  dif ferentiable  functions  of  x,t.  Part  (a) 
follows  from   these   two    observations. 

To  prove  part  (b) ,  we  differentiate  with  respect  to  x,  and  then  t, 
the   relation 


-91- 

X  =  x+(e,t)    ,      e   =  B(x,t)    ; 
we   get 

1    =  3^x+   3^,0=  3^x+   3^e  +  3^xt    . 

Using  equation   (5.34)    to   express  3  ,^x+    as  V  we    obtain, 

3^B    =   -  3tx+/3j^x+   =   -  V  3^      . 

This      is      equation     (5.32),      equivalent      to      (5.27).        We   appeal  now   to 
Theorem  5.6  to   conclude   that    (b)   holds. 
Lemma  5.10.      For  n    >    0 


('•")  R»  (  i^  )    2  0 


/ 
0 


(5.38)  /         Im  (   I^  )    dy    >   0 


Proof:    By  (5.12),    Ref     ,    j]    is      nonzero     only      on      the      set     I.      By 

R(n  ) 

Corollary      5.4,      P(n )     alternates    sign  on   the   intervals    that   make  up  I, 

and  by    (5.8)    it    is   positive   in   the   rightmost   interval.        On     the      other 

hand,      by      (5.4)      and      (5.5),      R(ri )     has      the      same   property,    so   (5.37) 

follows. 

Condition  (5.20)  ensures  that  f   Im  f  .  , 1  du  is  nonzero  only   off 

•'O    ^R(w)'^ 

the   set  I.   By  Corollary  5.4,  we  see  that  P()j)  alternatively  crosses 
from  negative  to  positive  and  positive  to   negative   on   the   intervals 

(32k+l,32k)  ^"^     such  a  way  that  it  is  positive  on  (B^,™).   Since,  by 
(5.4)  and  (5.5),  iR(y)  alternates   sign  on   those   Intervals   and   is 
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positive  on  (Bj,<»),  we  conclude  that  Im  (  .^  J  crosses  from  positive  to 
negative  on  the  Intervals  (3  2k+l'^2k^  ^"^  ^^  positive  on  (3  ^  ,«> ) .  From 
this   fact   (5.38)   follows,    completing   the   proof   of    the   lemma. 

To  construct  the  solution  of  the  variational  problem  (2.16-2.18), 
i|)*(x,t),  we  integrate  i|'*(x,t)  as  given  by  (5.16).  Theorem  3.13 
requires  that  \i)  (x,t)  =  0  where  nx  -  4n-*t  -  e^(n)  >  0  ;  in  particular, 
it  follows  that,  for  fixed  n  and  t,  i|;*(x,t)  =  0  for  x  large  enough. 
Thus  we    set 


(5.39)     .|;*(x,t)    =  -/   >l.*(y,t)    dy  =   /   Re[I^)dy, 


Lemma  5.8  shows  that,  for  fixed  n  and  t,  the  integrand  is  zero  for  x 
large.  Lemma  5.9  shows  that  d  ^ili  (x,t)  as  calculated  from  (5.39)  equals 
i|;*(x,t)  as  given  by  (5.17).  Inequality  (5.37)  of  Lemma  5.10  shows  that 
\li  is  a  positive  and  decreasing  function  of  x.  In  the  x,t  region  where 
X   <  x_(n,t).    Lemma   5.8  shows    that 

x+dn  ,t) 
x_(ri,t) 


and 


3xt*      =     ^t**      =      0   • 


Since  we   have  shown  i|;    (x,0)    =  (j)    for   x  <  x_(n  )  we    conclude 


x+(n ,t) 


P(n) 


(5.40)  /    ^fM-dy  =  <t,(n). 

x_(n  ,t) 


JL 

Finally,  since  we  have  shown  above  that  \|;  (x,t)  given  by  (5.39)  is 
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a    positive,    decreasing   function  of   x,    it  follows   from   (5.40)    that      0     _< 
\|)      <^  (J) ,    i.e.      that  \\)      belongs    to   the  admissible   set  A. 

The  main   result   of    this   section   is 

Theorem  5.11.      The  solution  of  the   variational   problem   (2.16-2.18) 
is   given  by 


(5.41)      ii;*(x,t)      ={ 


0 
x+(n  ,t) 


for 


P(n) 


X  >  x+(n,t) 
/      ^^ '^y       ^°'"  ^-(n,t)  <  X  <  X4.(n,t) 
<t>(n)  for  X  _<  x_(ti  ,t) 


Proof:    We     will     show      that      i|;    (x,t)      given     above      satisfies      the 
variational  conditions    (3.34).      Introducing 

(5.42)  e^(n,t)  =n  x+di.t)  -  L^)  (x^.(n  ,t)  ,t) , 

we   integrate   (5.18)  with   respect   to  x  to  obtain 


(5.43)      L^*(x,t)=  • 


,P(U)- 


nx-e+(n,t)    -         j  j    Im(-^)dn    dy   for  x>x+(n  ,t) 

x+(n.t)     0        ^(^'^ 


nx-e  ,(n,t) 


x_(n,t) 


for   x_(ri  ,t)  _<  X  _<  x+(ri  ,t) 


,p(m)- 


nx-e.(n,t)+  /  /    Ira(_iii4)dM    dy      for  x  <  x_(n,t), 

0  ^^""^ 


Differentiating   (5.42)    and  using    (3.50)    shows 


(5.44) 


•k  it 

8^e+  =  n3tx+-  L4)x(x+,t)    3tx+-  Li|;t(x^,t) 


=  4n^  . 
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We   know  from   (4.16)    that 


e+(n,0)  =n  x+(n  ,0)  -  l^   (x^(ti,0),0) 

=  nx+(n)  -  U)   (x+(n),0)  =e+(n). 


Therefore  integrating  (5.44)  we  see  that 


(5.45) 


e.(n,t)  =  4n^t  +  6,(n) 


and,  using  the  definition  (1.23)  of  a,  that 


(5.46) 


a(n  ,x,t)  =  nx  -  e^(n  ,t), 


Thus,  using  (5.46),  fornula  (5.43)  becomes 


(5.47)  a  -  Lij;*= 


X     n 


,P(li), 


/     /  Im(-^)dn  dy  for  x>x+(n,t) 
x+(r,,t)   0    ^(^^^ 


0 
x_(n,t) 


for  x_(ri,t)  £  X  £  X 


.P(vi). 


-/     /  Im(^il4)dy  dy   for  x  <  x_(n  ,t) 
0    ^^ 


Using  (5.38)  of  Lemma  5.10  we  conclude  from  this  that  ^  (x,t)  given  by 
(5.41)  satisfies  the  variational  conditions  (3.34).  This  proves 
Theorem  5.11. 

We  may  obtain  explicit  formulas  for  u(x,t),  u^(x,t)  and  u-^;(x,t) 
by  substituting  ij;^  given  by  (5.16)  and  ip^  given  by  (5.17)  into  (3.59) 
of  Theorem  3. 17: 
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(5.48)  u  =1  (n,^*)  =  --/  n  J^dn 


(3.3»  ;:^=-^(„3.,.,    ..^;       3QU. 


dn 


These  integrals  can  be  evaluated  by  contour  integration.  Silt  the 
plane  along  the  intervals  of  the  set  T  given  by  (5.1);  the  integral 
along  each  interval  is  then  equal  to  half  the  integral  of  a  clockwise 
contour  around  its  slit.  By  Cauchy's  Integral  theorem  we  can  deform 
the  contour  to  a  large  clockwise  circle  around  the  origin  containing  I. 
Then  letting  the  radius  of  this  circle  go  to  infinity  we  can  evaluate 
the   integral   by  computing   its    residue  at    infinity. 

For  large    \c,  |  we    have   by   (5.4) 

n 

r2(c)  =  1  I  (ej  -  c^) 

k=l 

n  ^2 

(5.52)  =    (-  c^)"    I     1(1-^) 

k  ^  ij  ^ijk  C  1^1 

where   the  sums  over  more  than  one  index  are  understood  never  to  repeat 
an  index  value  in  any  term.   Since 


(Hto)-l''2  =  1-|^+|^2_^^3+  0(0,^), 
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we    see   from   (5.52)   and    (5.5)    that 


^  =  i.n4(p&-LM|(Ia^^-i(^I.|e?))-L 


(5.53) 


^'^    k  ^    k        ij  ^  ijk  ?^        kl 


From   (5.8)   and   (5.9)  we   see    that 


P(o  =  ell  +4  +  4+  o(-4>]  . 


c2     c"^         \a^ 


Q(C)      =     cn-H2[,2-6(    l3^]-L  +  I|  +  I^H-0(-^)]. 

k         c^      c^      c^  ici^ 


Combining  this  with  (5.53)  gives 


^-l-("^4(p^)^ 


(5.54) 


k  k  i  j        C       IC  I 


g|j-l[l2,2+(V2  4(pa^-6(j6?6?))-L 
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(5.55)  +(y4+t(    I   3^)   Y2+j(    I   4)^ 

k  k 


k  ij  ijk  c^  Kr 


Using      (5.54)      and      (5.55)      to      find      the     residues,      we      compute      from 
(5.48-5.51): 

(5.56)  u   =1^  ^  I|^d,      =      -il   e2)_2a2 

(5.57)  J=   -^P,3  P(0  ,^   __  ^    ^   32^2   .4   ^    ^   3?3?)   -^I   3^)«2-f«4 

k  i  j  k 

(5.58)  7.-J-f,,«gd,      -      (I6j)2-4(l6?6?)+i.2 

k  ij 

(5.59)  ^  =  i^/)c3^dc   =   -[    I   e^)^  +  4(    Ul)    [I   ,1,]] 

k  i  j 

ijk  ■"       k 


Note    that   together,    (5.57)   and    (5.58)  yield  an  algebraic   identity 


T2   =   M    ^h]  +   ^[    I   &lh2  +  ^4 
ij  k 


Now  we  consider  the  region  of  x  and  t  where  n,  the  cardinality  of 
B(x,t)  is  equal  to  one.   In  that  case  formulas  (5.56  -  5.59)  become 

(5.60)       u   =  -  e^  ,   u2  =  e^  ^   ^3   =  _  36  ^ 
which  are  (4.39)  and  (4.41). 
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6.    Asymptotic  Behavior   for  Large  Time 

The  main   result   of   this   section   concerns   the   values    of  u(x,t)   as    t 
-y   00  : 

Theorem  6.1.      Let  6    be  any   positive   quantity: 

(a)   For   x/t    <  -  6 


(6.1) 


u(x,t)    =  0(t"2), 


(b)    For  6    <   x/5    <  4-6 


(6.2)  u(x,t)    =  - 


^        (a(x,t))   +o(t"b    , 


2TTt 


where 


(6.2)' 


o(x,t)    =2.  (x/t)l/2    . 


(c)   For  x/t   >  4  +  6 


(6.3) 


u(x,t)    =  0(t    ^), 


The      proof      is    based  on   the    variational  characterization   of  u .      We 
use   condition   (3.35): 


(6.4) 


ijj*(n,x,t)   = 


0  where      nx-4Ti  ■^t-e^(Ti  )    >   0 


(n)     where     nx-4n ^t-e_(n )   <  0. 
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According  to  Lemma  1.4,  e^(n)  and  S.Cn)  are  continuous   functions   that 
vanish  at  n  =  0.  Prom  this  and  (6.4)  we  deduce 
Lemma  6.2.   (a)  For  x/t  <  -6 , 


(6.5) 


i^*(n)  =4>(n)  for  O(t-l)  <  n  <  1  . 


(b)  For  6  <  x/t  <  4-6 


(6.6) 


li*  (n)  =  { 


0   for  0(t"^)  <n  <a(x,t)  -0(t"l) 


(n)  for  a(x,t)  +  0(t  ^)  <  n  <  1. 


(c)  For  x/t  >  4-W 


(6.7) 


i|;*(n)  =  0  for  0(t"^)  <  n  i  1. 


We  show  now  that  parts  (a)  and  (c)  of  Theorem  6.1  follow  from 
parts  (a)  and  (c)  of  Lemma  6.2.  It  follows  from  (6.5)  and  (6.7)  that 
In  both  cases  the  parameters  g  .  Introduced  in  Section  4  lie  in  the 
range  (0,  0(t~  )).   We  assume  that  there  is  in  fact  only  a  single  ^. 

According  to  formula  (5.60) 


u  =  -  e^ 


since   3^  _<  0(t~^),    (6.1)   and    (6.3)    follow. 

We   turn  now  to   the   interesting  case      (b);      here      it      follows      from 
(6.6)    that    the  parameters  ^.   are   contained   in   the   intervals 
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[O,0(t"b]  or   [a(x,t)  -O(t-l),  a(x,t)  +0(t"^)] 

with  at  least  two  3 's  in  the  second  interval.   We  assume  that   in   fact 
there  are  exactly  two,  3^  and  ^2  >  i"  the  second  interval, 
and  one,  3 ^  ,  in  the  first.   Thus  (6.6)  can  be  rewritten  as 


(6.8) 


Here 
(6.9) 
and 
(6.9) 


0        for     n      in   [33,32! 
^*(ri  )    =     i    'j)(n)     for     n      in   [3i,l] 

in   (0,(J)(n))  for  n    in  (0,33)   and    (32. Bp. 


13  j   -  o(x,t)|    <  0(t    b,      j    =    1,2 


3  =  o(t  h. 


Ihe  next  lemma  nails  down  more  precisely  the  distance  between  3i 
and  32'  We  assume  that  <^  is  continuous  in  (0,1)  and  that  6^  is 
dif f erentiable . 

LemnH  6.3. 


(6.10) 


log(3i-32) 
lim  f— (})(o(x,t))  =  -  2ir 


t+«> 


Proof:    According   to   the   variational   condition   (3.34), 


(6.11) 


Lij;      =  a 
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at    every  point  r\    where  0   <  i|j    (n  )    <  <t>(n).      In    particular,      according      to 
(6.8),      (6.11)   holds    for  n    in  (32. 3i).    and   so   by    continuity  also  at  n    = 
&^    ,    i  =    1,2.      We    now   set  n    =  3^^   in   (6.11);    using   the   definition   (2.20) 
of  L  we    get 


2tt  i     ^  ^e^-Hj- 


(6.11)'     a(6.)    =    (L,j;*)(6p    =  J-/    log  (^_)2^*(u)    dp    ,       i  =    1,2. 


According   to   (6.8),   ijj    (y  )    =  0  for  3  3  _<  y  _<  g  2  ^""^  ^    ^^"^    "^  <^(v)   for  g^  _< 
y    <   1.    Using   this    information,    we   obtain   the   following: 


(6.12)  a(ep    -a(B2)    =    (L'l'    X^i)    "    (L'l^    )(B2) 


1   ^3  (ei-P)(6  2^)  2  *,  ,  , 

2tt  'q  (e2-u)(ei-Hj) 


1       ^  6?+M     o       * 


'2 

1 


+  1,1    log(^)2/(p)    dp 


^2 

1 


+  ^/    log  [^^^)'^  *(u)    dp 


2^  e^"'°  ^e2ni- 

We    claim   that    the   first    three    terms    on   the   right   are  0{&^-Q2)- 

This    is   obvious    for   the    first    two   terms;    in   the    third    term  we   introduce 

as   new  variable   of   integration 


V  -  e>2 

X  = 

^1-  ^2 
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We  get 

1 


(^1-^2)  /  log  (^-T^)  1^*  '^ 


Since  ij)   <(}>_<  const.,  the  above  quantity  is  indeed  0(3^  -  62^*    Thus 
we  can  rewrite  (6.12)as 


1  1      ^I'V'  9 
(6.13)     a(e,)  -  3(30)  =  0-  /  log  L )  't'(M)  dp  +  o(e,-e2) 

^1 


We  further  split  the  integral  on  the  right  as 

/    =     /      +   /    . 

Y   a  small  quantity.   For  fixed  y,    the  second  integral  is  0(3,^-32)  '  1" 
the  first  integral  we  write 

4>(p)   =  (l)(Bi)  +  e(y)  ; 

e  (y )   is    small   in   [B^jg^+y]   when  y    is   small.      We    Introduce 

<   =  M    -  S>i 

as    new  variable   of    integration;    this   allows   us      to      rewrite      the      right 
side   of   (6.13)    as 


(6.13)'  ij    [<t,(B    )4^]    iog_^-l__  d<  +  0(B^^2) 

TT      ^  KT(,P  j^-g  2.) 

=  1   [4,(3  P+e]    (Bi^2>    log   (B1-B2)   +0(Pi-B2) 
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We      turn     now      to   the    left   side   of   (6.13)    and  apply    the  mean  value 
theorem;   we    get 


(6.14)  a(ep   -3(62)    =\^^l-^2^ 


where 

(6.15)  ^  =  a^(e"),      ^2   <  e"  <  3i 

It  follows  from  (6.9)  that 


(6.16)         3  =  a(x,t)  +  0(t  M 


We  use  now  definition  (1.23)  of  a(n,x,t)  to  write 

a  =  X  -  12tn^-de^./dn  . 


using  this  and  (6.16)  in  (6.15)  we  get 

7^  =   x-12tx/4t-K)(l)  =  -2x  +  0(1)  . 

Using  this  in  (6.14)  we  get 

a(6i)-a(e2)  =  -2x(B^-02)  +0(31-^2) 
This  is  equal  to  (6.13)';  dividing  by  (3^-02)  ^e  get 


-2x-K)(l)  =  i  ((t>(0i)4t)  log  (Bi-62)  +0(1) 
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Using  (6.9)  once  more  we  note  that 

♦  (Si)   =  <t)(a(x,t))  +  e(t)  , 

where  e(t)  -»■    0  as    t  ->■  <»    ;    relation   (6.10)    of  Lemma   6.3   follows. 
We    make   use   of    fortiula   (5.56)    for   the    limit    solution: 


(6.17)  u   =   -{^l-^l-^])    -  2a2    ; 


02   is   a   function  of  6]^,B2»B3  defined   through    (5.20).      With  n  =    3, 


P(z)         =       Z-^     +   023 


so   (5.20)    is 


^2  n  ^  +  ct  9n 

/     ,    r     dn  =  0 

R         R(n ) 
33 


Thus 


(6.18)  a,    =   -  j(3^    /  j(^>     . 


where 


(6.19) 


and 


''"■trnj'"-   '■'•'• 


(6.20)  R2(n)  =  (n^-e|)(e|-n^)(B5-n2) 


According   to   (6.9),    (6.9)',    as    t  +   0,   61-B2  ■*"    '^  ^nd     B3     +      0.      We 
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shall        now     determine      the      asymptotic      behavior      of     02     under      these 
circumstances . 

Lemma    6.4.      As  3^-62  arid  63  >    0, 

(a) 

(6.21)  efj^^^  -  J^^^  -  Bi  *  0 


(b) 


(6.21)'  j(l)    ~     J-   Hog    (Bi-B2)l 


Proof:    (a)   Using    (6.19)    and    (6.20)    we   get 


(6.22)  b2j(1)_  j(3)=  j      4r^d^ 

1       R(n ) 


:2 


,2__2 


^1"^  ,1/: 


=  /     '1 -{-^ ]''^  dn 

B3   (r^^-Q2^  ^2^^ 


The  integrand  on  the  right  in  (6.22)  tends   to   1   uniformly   on  every 

closed   subset   of   (BtjPo^*  ^^'^      ^^      dominaed  over   that  interval  by 

1  /? 
const ./ [  (n-P  3)  (B  2"^  )  ]     »  ^^   integrable   function.    Therfore   by   the 

principle   of   dominated   convergence   as   B  3  ->•  0  and  3  o"^  1  "*"  ^»  (6. 22) 

approximates  B 1  ,  as  asserted  in  (6.21), 

To  prove  (6.21)'  we  separate  those  factors  of  R  which  are  small  at 

the  end  point  B 1  : 
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The  value  of  the  first  factor  at  n  =  Bo  is 


23?  +  0(3i  -  B2)  +  0(^3) 


From  this  we  deduce  easily  that 


(6.22)'  jO)-  J_  /^  ^ ; 

^^1  B3  [(ei-n)(3  2-n)]^^^ 


elementary   calculus  shows  that  the  integral  on  the  right  in  (6.22)'  is 
■~  -logCgj-g^)  ;  this  completes  the  proof  of  part  (b)  of  Lemma  6.4. 
Using  (6.21),  (6.21)'  in  (6.18)  gives 


r2j(1)_j(3) 
ao  =   -  J^3)/j(1)   =  ^1^ 
"2  ^'  -TT) 


26  2 


1   log(6^^2)    ^log(Bj-62)^ 


Setting   this    into   (6.17)    and   using    (6.10)    we   get 


26 

(6.23)  u(x,t)    =  6?   -  60    -  6?   -  — -  (t)(o(x,t))   +  o(t~b. 


It    follows    from   (6.20)    that 
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Using  this  to  estimate  the  first  two  terms,  (6.9)'  to  estimate  the 
third  terra  and  (6.9)  to  estimate  the  fourth  term  we  obtain  (6.2).  This 
completes    the    proof   of    Theorem  6.1. 

We  conclude   with  a  number  of   observations: 

(i)   A   similar  argument,    baed   on    formula    (5.57)    for  u    ,    yields 

(6.24)  u2(x,t)    =  _^<t,(a(x,t))   +o(t"b. 

6Trt^ 

The    first    relation   of   (4.42)is 

(6.25)  u^   =    3u2^   ; 

the  leading  terms  in  the  asymptotic  descriptions  (6.2)  and  (6.24)  are 
consistent  with   this    relation. 

(il)      Relation   (6.6)   suggests    the    following  crude  approximation   to 
^    (n,x,t)    in   case   (b)  : 

\|;*(n,x,t)    =  H(n    -  a(x,t))  c})  (n )    , 
where 

H(vj)    =0     for  vi    <   0   ,      =1  for  p   >   0     . 

Then,    using    (6.2)', 


*  1 

^Ar]  ,x,t)   =  -  6(n   -  a(x,t))  __  ,j)  (n  ) 

4(xt)^/2 
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Using    this    expression   in   formula    (3.59)    for  u  gives 


""/      ^N  ^    f      ,*\  ^      /      ^N   <t>  (a(x,t)) 

u(x,t)      =-   (n,t„)    =  --0(x,t)  1 '    ■„ 

4(xt)l/2 


<t>(a(x,t))    , 


2iit 


in  agreement   with   the   leading   term  in   (6.2). 

(iii)  Next  we  show  that  (6.2)  can  be  derived  from  the  zero 
dispersion  limit  of  Tanaka's  description  of  the  large  t  behavior  of 
solutions  of  KdV.  Tanaka  shows  that  for  large  t,  and  for  0  <  x  <  4t 
every  solution  is  a  wave  train,  i.e.  approximately  a  superposition  of 
solitons 

N 
(6-26)  u(x,t;e)  >   I    six  -   UT\lt  -  &^    ,  n^)    ; 

1 

here 

(6.26)'  s(x,n)    =     -   2ri      sech   (nx/e)    . 

The   width   of    such   a   solition   is 


(6.27)  const.   £ /n 


It    follows    from  the   precise   form   (1.18)^^  of   Weyl's    law   that 


(6.28)  n^+i    -^n     -     -^^^  '      ^n     ^^    f^n'^n+l^       • 


Peaks  of  the  wave  train  (6.26)  are  located  at 
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(6.29)  Anjt  +  6^   . 


Using    (6.28)   we   see   that    for   t   large  the  peaks      are      separated     by     the 
distance 


(6.29)'  '  .      ai^t(n„+i   -nj   +6^+1   -6^  - 


Sirrij^tE 


Comparing      this      to      (6.27)    we    conclude  that    for   large   t    the   individual 
solitons    in   the   wave   train  are   well  separated. 

The  wave   number  n    of    the   soliton   that   peaks   at  x  at   time   t      is      by 
(6.29)   and   (6.2)', 


^n     _ 


-     o(x,t)    =      (x/4t)l/^ 


provided      that      t      is      large  and  0<  x/t    <  4.    Setting   this    in    (6.29)'  we 
conclude   that   for   t    large   the    density   of   solitions   at    x  is 

(6.30)  _     (J)(a(x,t))/4iTe(xt)^/2    ^ 

It    follows    from   (6.26)'   that    the   area  under  a   soliton   is   4ne    ;      at      the 
wave   number   (6.29)    this   is 


(6.30)'  2(x/t)^/2  e    . 


The   asymptotic  area  density  in  the  wave  tain  (6.26)  is  the  product  of 
(6.30)  and  (6.30)': 
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(2TTt)  ^  (})(a(x,t))  . 

Since  the  asymptotic  area  density  is  the  weak  limit  u,   we   obtain  yet 
another  derivation  of  (6.2). 

(iv)  In  deriving  (6.2)  we  have  assumed  that  the  function  (^  is 
continuous.  When  the  potential  has  several  local  minima,  the  function 
(})  has  discontinuities;  in  these  cases  Theorem  6.1  remains  true  as  long 
as  a(x,t)  is  bounded  away  from  the  points  of  discontinuity  of  ({) . 
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7.    Initial   data   with  u(-°°)    ^  u(<»). 

In    this   section  we    study    solutions   of 

(7.1)  u^   -    6uu^+e2u^^^  =    0 

whose   initial  value  u   satisfies 


(7.2)  lim   u(x)    =   -1,      lira  u(x)    =   0, 


and  is  an  increasing  function  of  x. 

The  solution  of  such  an  initial  value  problem  can  be  obtained  as 
the  limit  of  solutions  whose  initial  data  u  satisfy  the  conditions 
laid  down  in  Section  1; 


(7.3)  lim  u^(x)  =  u(x)  , 


n-H» 

uniformly   for  x  >  Xq    ,      xq      any     number.        Denoting      by     Uj^(x,t,e)      the 
solution   of   (7.1)   with   initial  data   u   (x) ,    it  can   be    shown   that 


(7.3)'  lira    Uj^(x,t;e)    =u(x,t;e) 

n-Hxj 


Although  for  the  solutions  of  KdV  equations  signals  propagate  with 
infinite  speed,  the  dependence  of  the  solution  at  fixed  x,t  on  the 
values  of  the  initial  data  at  y  diminishes  to  zero  as  y  tends  to  », 
uniformly  in  e  .  Therefore  we  conjecture  tha  the  limits  n  ->•  «>  and  e  -»■  0 
can  be  interchanged: 
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Conjecture   7.1: 


(7.4)  d-lim     u(x,t,e)    =lim   d-lim     Uj^(x,t,e). 

£+0  n-Hx>        e  +  0 


The      Inner   limit    on   Che    right    in   (7.5)    is    characterized    in  Theorem 
2.10  in   terms    of    the   solution   of   a  minimum   problem    (2.16):  : 


(7.5)  Uj^  =  d-lim  u^(x,t,e)    =  8^Q^    , 


where 


(7.5)'             Q*    =  min      Q^(^;x,t), 

and 

(7.6)         )    q    (xl,)  =-   (a^,rp)    -  -   (L^.ij;) 

"  TT             "                      IT 


By   (1.23) 

(7.7)  a^  =  xn  -  4kn3  -  ei")(n)   ; 

by    (1.25) 

xi")(n) 

(7.8)  (|)^(n)   =        J         1 ^  dy    , 

x(n)(n)    l%(yHn2|l/2 

and  by   (1.  16) 

00 

(7.9)  e(")(n)   =  nxi")(n)  +    j      (n   -  (u„(y)  +  n^]  ^^^  ^y 

x(n) 

xj^)(n)   are   defined   by   ( 1 .      )  : 
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(7.10)  Uj^(xi"^<^'^))    =  "  '^ 


Since   by   (7.3),    u^  tend   to   u,    and    u  satisfies    {1.1),    it   follows   that 
(7.11)  lim  xi")(n)  =  -0°  ,     lira  4"''(n)   =  x^.(n) 

From  this   we    deduce   easily    that 

•  lira  ei'^)(n)  =  e^(n)  , 
n->«> 

lira  aj^(x,t,n)    =a(x,t,ri)    , 
n-HJo 

lira  ())j^(n )    =  "    ,  ■  ' 

n«% 


uniformly    in  n  •      It    follows    from  this    that    for   fixed  i|;    in  L    , 


(7.11)  lira  Q^(,)j;x,t)    =  Q(;|;,x,t) 

We    surmise    tha   the   limit    of    the  minima  Q      defined   in  (7.5)      can  be 
characterized   as    the   minimum   of   (7.11): 
Conjecture  7.2: 


(7.12)  .  lim  Q*   =   Q* 


where 


(7.13)  '  Q*   =  rain  Q(,j;)    , 

(7.13)'  Q(ij;)    =  -   (a, I);)    -  -   {h^  ,Mf) 

TT  11 


Combining   conjectures  7.1   and  7.2  we    deduce 
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Conjecture  7.3: 


(7.14)  d-lim       u(x,t,e)    =  d^Q*    , 

where  Q      is    defined  by   (7.13). 

It  is  far  from  obvious  that  the  minimum  problem  (7.13)  has  a 
solution.  Nevertheless  it  follows  from  the  positive  definiteness  of  Q 
that  it  ^  satisfies  the  variational  conditions  for  the  minimum 
problem,  then  Q(i|^)  achieves  its  minimum  at  i|j  ,  and  only  at  ^  .  The 
variational   conditions    (3.34)    in  the    present    case  are: 


(7.15)  ^li*   =   0     where      a  -  L^^)*   >   0 


(7.15)  '  a   -  LiJ;*   >  0     for  all  n    . 


We   turn  now  to   the    special    initial    function 

[    -1      for     X   <  0 

(7.16)  u(x.)    =        ] 

[    0        for     X  >  0 

In   his    case    x^(n )  =    0,    6+(ri)   =    0,    and 

(7.17)  a   =  xn  -   4tn^    . 

Equation  (7.13)    can   be   solved   explicitly: 

Theorem  7.4:      For      a     given     by      (7.17),      the      solution     ij;      of    the 
variational    problem   (7.13)    is  as    follows: 

(a)   For 

(7.18)  X   <  -  6t  .  ■  .  - 


(7.18)'  t*   =   ^^^^-(^<^-^)^ 

(l-n2)l/2 


(b)    For 


(7.19)  -  6t   <  X    <   4t 


0  for  n   < 
(7.19)'  ^*  = 


^     for  B   <  n    <  1 
^0 


(7.20)  R2(n)    =   (i-n2)(^2_g2)    ^ 

(7.21)  S(n,x,t)    =   12tQ   -  xP   , 

(7.21)'  q(ti)   =  n^  +  Y^n^  +  Y2   ,     P(n)=n^+a 
where 

(7.22)  a    =  -  l(2)/i(0) 

(7.23)  Yi  2 


1    +  6^ 


(7.24)  Y2    =   -(Yil^2^    +1^^^    /   1^°^    . 

I*-^-'   denotes    the   complete   elliptic    integal 

(7.25)'  R2(n)    =    (l-n2)(B2-^2)    ^ 

The    parameter  3    is    related   to   x/t    through 
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(7.26)  S(B,x,t)    =   0  . 

(c)   For 

(7.27)  4t   <   X    , 
(7.27)  '                                ,];*  =    0 


Proof:  Vie  verify  (c)  first.  Clearly  \\)  is  nonnegative  and  so 
admissible.  We  show  now  ttiat  the  variational  conditions  are  satisfied. 
Since  ij;      =   0,    hli      =0;   combining   this  with   (7.17)   we    obtain 


(7.28)  a    -  Lij;*     =      xn    -    4tn  ^    . 


Clearly,    if   (7.27)   holds,    (7.28)    is  >    0   in   the    range   0   <_  n     _<      1;      thus 
(7.15)    and    (7.15)'   are   satisfied. 

Next      we      take      case      (a).      It    follows   from    (5.16),    (5.17),    n  =   1, 
that    for  ^^^*   defined  by    (7.18)', 


a   -  Lijj*   =   0     for     0  <  n    <   1. 


This  shows  that  the  variational  conditions  are  satisfied.  It  further 
follows  from  (7.18)  that  i^*  as  defined  by  (7.18)'  is  >  0  for  n  >.  0  and 
therefore    is   admissible. 

We  turn  now  to  case    (b) .      Fix  g    in   [0,1],    and   define      x/t      through 
the    relation  (7.26). 

"k  '  -  ' 

We  can  write   \\t    ,    given  by   (7.19)',    as 


(7.29)  i|;*    =  »(/t    -    12tiJ;^ 


-117- 


where 


^1 


(7.30) 


,  * 


0  ■      ■  for  n    <  B    , 

—  ,      P  =  Ti^  +  a        for  B   <  n    <  1 
^0 


0      ,  ■        ■  -  for  n    <  B    , 

-^  ,     Q  =  n^  +  Y^n^  +  Y2      for  B  <  n   <  1   , 
^0 


see    fornulas   (5.16),    (5.17),    n  =    2. 

We    see   that  ij; ,    and  ^2  ^^^   real   parts      on      the      real      axis      of      the 


ana 


lytic    functions,    see   (5.16),    (5.17),    n   =  2; 


f    ,  s        P(n) 
^  Ro(n) 

(7.31) 

f   /  N   _   Q(n)        .2 

Clearly,      fj^      is      bounded      in  the    upper   half    plane;    using    (7.23)    we   see 
so   is   f2. 

Using    (7.29),  (7.31)   and    (7.21)    we  get 


that    so   is   f2. 


(7.32)  ^*   =  Re    I    , 

(7.32)'  •         f(n)    =  1^+  i(x-12tn2)    . 

Ro('^) 


Using  relation   (3.41)   between   the    operator  L  and   the   Hilbert      transform 
H   we    deduce    from    (7.32),    (7.32)'    that 
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(7.33) 


L\()*   =  J    H  4)*   dvi 


Im(f)    dy 


-  j    i^du   +xr, 


RTTT 


-  4n^  ,     n  < 


-  J  1^  dp  +  xn  -   4n^  ,     n  ^ 


0 


RTTy 


where  R    ,   defined  by   (7.25)',    is   positive   on  (0,B). 
Relations   (7.22)    and   (7.24)    guarantee   that 


(7.34) 


Combining   this  with   (7.33)    and  using    (7.17)    we    obtain 


(7.35) 


L\l)' 


i    I-  dp   +  a 
0  ^ 


for     e    < 


for     n    > 


We  are  now  ready  to  verify  the  variational  conditions.  Since  R(ij) 
is  of  one  sign  in  (0,0),  (7.34)  implies  that  S(p)  has  a  zero  in  (O.g). 
According  to  (7.26),  S(3 )  =  0;  since  S  is  a  quadratic  polynomial  in  u^, 
it  can  have  no  more  than  these  two  roots  on  the  positive  axis  and 
furthermore  these  two  roots  are  simple;  we  draw  now  some  consequences: 

Since  the  leading  term  in  S  is  12tp'^,  S(p )  >  0  for  y  large.  It 
follows  therefore  that  S(p)  >  0  for  u  >  3;  then,  since  S  has  a  simple 
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root  at   3.    S(m)    <  0  for  e    <  y   and  near  3.      Then,    since     S     has     exactly 
one   simple    root    in  (0,3),    S(u )    >   0  for  p    small,    see  Fig.    2: 


Figure  2 


Since  R   is  positive  in  (0,3),  it  follows  that  for  n  small  and  > 


0, 


-120- 


(7.36)  J  l^i-)  dy  >  0  . 

0  ^^^^ 


We  claim  that  (7.36)  holds  for  all  n  <  3  ;  for  clearly,  it  holds  for  n 
<  Hq  ,  where  n q  is  the  root  of  S  in  (0,g).  Since  S  is  negative  on 
(Hq.B),  we  deduce  from  (7.34)  that  for  n  between  Hq  and  3, 


J  I  dM  >  j  I  dy  =  0 
0  '^      0  '^ 


This    proves    our  contention. 

Combining    (733),    (7.3A)   and    (7.36)    we   have 

>  0      for     n    < 

T  ,  * 
a   -    Li|i 

0     for     B    <  n 

In  view  of  (7.19)'  this  shows  that  the  variational  conditions  (7.15), 
(7.15)'  are  satisfied. 

We  note   finally  that,  since  S(p)  >  0  for  p  >  g ,  It  follows  from 

(7.19)'  that  ij;   >  0  and  therefore  admissible. 

To  complete  the  proof  of  part  (b)  ,  we   have   to   show   tliat   as   (3 

varies  from  0  to  1,  x/t  varies  between  -6  and  4.  From  (7.21)  and  (7.26) 
we  can  express 


(7.37)  _  =   12 
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X        .„  e^^  +  ^l^^  +^2 


t.  b2+^ 


The  coefficients,  y^  andY2  ^^^^  given  by  formulas  (7.22),  (7.23)  and 
(7.24)  in  terras  of  the  complete  elliptic  Integrals  I^-J-'  defined  in 
(7.25).      By   elementary   calculus  we   deduce 

Lemma    7.5:    As   3   ^   0 

(i)        l(0)(3)=^/2 

(ii)      l(2)(B)4p2 

4  ,  -  - 

(ill)  i(M(e).^34 

16  ,  ,.,-■,  ■ 

We   deduce   from   this    fact    that    as    6   ->•   0 


a     ^ 


B^  _   1  g2 

X  '       ^1    ""         2"  '      ^2    ==    — 


Setting    this    into    (7.37)   gives 

(7.38)  21   =     12  iL/i  -     -  6     as  6   ^   ) 

•^  b2/2 

Using   intermediate   calculus   we    deduce 
Lemma    7.6:    As   S   ->•    1 


(i)        l^(&)^^   llogd   -  B)| 
(11)      I^O)(B)    -  I^2)(g)  ^1 
(ill)    l(2)(3)    _  i(4)(g)  ,1/3 
We    deduce    from    this    that 


-122- 


oc    =     -   1   +  Iq^    ,     Yi  =      -1    ,      Y2    =    ^  Iq 


3 


Setting    this    into    (7.37)   gives 


(l/3)Io^ 
(7.38)'  21    ,    12  —  =  A     as     B   >   1 

t  x-1 

^0 


Since  x/t  as  defined  by  (7.37)  is  a  continuous  function  of  g ,  it 
follows  that  as  3  varies  from  0  to  1,  x/t  varies  between  -6  and  4. 
Since  according  to  Theorem  3.12  for  each  x,t  there  can  be  at  most  one 
function  \|j  that  satisfies  the  variational  condition,  it  follows  that 
the  relation  of  g  to  x/t  is  one-to-one.  This  completes  the  proof  of 
Theorem  7.4. 

According    to   formulas    (5.56),    (5.57) 

u     =      _   1   _  g2   _   2a 


u2   =    (14^2)2   -  ^^2  +  4    (j^2)a 


Setting    the   definitions    (7.22)-(7. 24)    into    these    formulas   yields 

(7.39)  u     =  -    1   -  $2   +  2l(2)(e)/i(0)(g) 

and 


(7.40)  u2   =    1  +  e^   +  I.  ^2   _  A    (i+g2)x(2)(g)/i(0)(3) 
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Theorem   7.7:    (a)   For  x  <  -6t 

u(x,t)    =   1    ,      u^Cx.t)    =    1 

(b)  For     -      6t   <   X   <   4t,    u  and  u^   are   given  by    (7.39)   and    (7.40), 
where   g    is   a    function   of   x/t    defined  by    (7.37)   and    (7.22)-(7 .24) . 

(c)  for  4t   <   X 

u(x,t)    =  0    ,      u2(x,t)    =   0   . 
Note  that   weak    limits  satisfy 


(7.41) 


0   <  "?   <  u^ 


Using  formulas  (7.39),  (7.40)  we  get  amusing  inequalities  for  complete 
elliptic  integrals  from  (7.41). 
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